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ABSTRACT 


Possible improvements to the minimum resolvable temperature difference (MRTD) 
entered into TAWS are considered. FLIR92 is modified to include atmospheric turbu¬ 
lence that depends on height in the atmosphere. Resultant MRTDs are compared to 
the operational FLIR92 MRTD predictions excluding atmospheric turbulence. The 
difference in the MRTD is only apparent in the higher frequency regime and is less 
than 0.05% of the MRTD values for a typical test case. MRTD is calculated by 
FLIR92 and NVThermlP over desert and marine locations and the resultant MRTDs 
entered into TAWS to compare maximum detection range. NVThermlP yielded a 
larger maximum detection range by up to 1.5% over the desert and 2% over water. 
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I. 


INTRODUCTION 


Target Acquisition Weapons Software (TAWS) is used operationally to predict 
the performance of electro-optical weapons and their navigation systems. For sensors 
operating in the infrared (IR) spectral band, the minimum resolvable temperature dif¬ 
ference (MRTD) is determined from a desktop computer model such as the Forward 
Looking Infrared 92 (FLIR92) or the Night Vision Thermal and Image Processing 
model (NVThermlP). Both models are capable of providing an MRTD, which may 
be manipulated for use in TAWS, but there are differences in their calculations that 
become apparent at high and low frequencies. Also, although FLIR92 does not oper¬ 
ationally include the effects of atmospheric turbulence, NVThermlP does for limited 
conditions. 

In this study, FLIR92 and NVThermlP are compared in their current opera¬ 
tional forms for various atmospheric conditions. Then FLIR92 is modified to include 
atmospheric turbulence and the two version of FLIR92 are compared for the varied 
atmospheric conditions. 

All acronyms in this report are listed in alphabetical order in Appendix E. 

A. MOTIVATION FOR COMPARISON 
1. Current Status of TAWS 

Currently, the MRTD used in TAWS predictions of target detection by IR 
sensor is calculated by FLIR92 from basic system parameters. While this has been 
sufficient for previous iterations of TAWS, recent improvements in model resolution 
and the image resolution of sensors being modeled suggest that perhaps NVThermlP 
may be a better choice. NVThermlP incorporates a number of improvements intended 
to improve the determined MRTDs for the sensors in the extrema of low and high 
spatial frequencies, in other words for very large and very small targets. NVThermlP 
also includes optical turbulence, which should improve the MRTD calculation. With 
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increased sensitivities of sensors, including optical turbulence should improve the 
MRTD prediction. Since FLIR92 and NVThermlP behave differently at high and low 
frequencies in marine or desert environments, one model may substantially improve 
the MRTD calculation, and thus the TAWS predictions, over the other model. In 
particular, the temperature difference between the surface and the ambient air may 
influence model calculations and this temperature difference is much greater during 
the day over a desert. 

2. Differences Between FLIR92 and NVThermlP 

Both FLIR92 and NVThermlP are desktop computer models developed to 
predict standard summary performance hgures of merit for thermal imaging sys¬ 
tems. They are PC based programs that model passive sensors which detect emitted 
and reflected radiation. Using basic system-level parameters, these models calculate 
the modulation transfer function (MTF), the noise equivalent temperature difference 
(NETD), the MRTD, and the minimum detectable temperature difference (MDTD) 
for sensors looking at specific targets. FLIR92 and NVThermlP both predict the 
MRTD that a human can discriminate when using a thermal imager operationally. 
NVThermlP also predicts the range at which target acquisition will be successful, 
using the specihed thermal imager. Thus, the output from each model is used to 
determine whether or not a system will achieve the MTF, system noise, MRTD, and 
MDTD that is required to effectively perform a given mission. For each mission, the 
conditions necessary to meet the target acquisition and discrimination requirements 
may vary. 

FLIR92 and NVThermlP share certain basic assumptions. These assumptions 
are stated in the NVThermlP User’s Manual (NVESD 2001), and follow the steps of 
modeling FLIR systems outlined in the Infrared and Electro-Optical (IREO) Hand¬ 
book Vol.4 (1993). All MTFs are assumed to be separable, so the total system MTF 
is calculated as the product of all sub-system MTFs. This approach reduces calcula¬ 
tion complexities because the analysis is simplihed to one dimension and cross-terms 


2 



are eliminated. However, there is almost always some calculation error associated 
with assuming separability. This error becomes signihcant in certain cases such as for 
diamond-shaped detectors, but is neglected in all FLIR92 and NVThermlP calcula¬ 
tions regardless of detector shape and type (NVESD 2001). 

An MTF describes the spatial frequency response of a system. It is the contrast 
at a given spatial frequency compared to the contrast at a lower frequency. The total 
system MTF is a product of all component MTFs. Each MTF may be calculated 
as the Fourier transform of the point spread function (PSF) of the component, or 
the response of the imaging component to a point of light. An imaging system 
will experience some degradation of the image due to imperfections in the optics, 
electronics, or even the observer eye and this degradation is what the MTF describes. 
If some system component did not contribute to the degradation, then the MTF for 
that component would be unity (Driggers 1999). 

Both FLIR92 and NVThermlP assume all blurs to be symmetrical in order to 
keep all MTFs real. The models assume that there exists a region of the held of view 
(FOV) that is isoplanatic. Goodman (1968) explains that in an isoplanatic, or space- 
invariant, linear imaging system the image of a point-source will change location, 
but not functional form, as the point source moves. The approximately isoplanatic 
region of the FOV is modeled in FLIR92 and NVThermlP by a linear shift-invariant 
process and the MTF in that region is approximated by symmetrical blur. The blur 
is the result of real world effects on the image and can be due to aberrations and 
other manufacturing defects in the optical system. For a point source, the image 
is called a blur circle because diffraction, aberrations, manufacturing defects, and 
assembly tolerances prevent perfect resolution of the singularity (Driggers 1999). In 
FLIR92, the blur MTF is geometric and is discussed in further detail in the next 
chapter. The symmetrical blur approximation is not an accurate rehection of real 
world systems, so the blur approximation for the optics is not completely correct. 
Similarly, assuming symmetrical blurs does not accurately rehect the electronics used 
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in the various sensor systems under consideration: a low pass filter would not result 
in a time delay or phase shift using this assumption, for example (NVESD 2001). 

NVThermlP provides some new capabilities beyond FLIR92. NVThermlP 
provides target acquisition performance predictions for staring imagers, not just the 
Erst and second generation thermal scanning sensors. In the NVThermlP calculation, 
the MTF representing the function of the human eye includes factors that are ignored 
in FLIR92. The MRTD prediction is changed from FLIR92; NVThermlP produces 
contrast transfer functions (CTFs) as the primary output. A laboratory MRTD is cal¬ 
culated, but is not comparable to the MRTD output from FLIR92. Also, NVThermlP 
now uses the Targeting Task Performance (TTP) metric to predict the probability 
of target detection, recognition, and identification at given ranges. FLIR92 uses the 
Johnson criteria, but does not predict probabilities of detection, recognition, or iden- 
tihcation (NVESD 2001). In this report, the TTP capabilities in NVThermlP are 
neglected since only the MRTD results from FLIR92 are compared to those from 
NVThermlP. 

3. Advantages and Disadvantages of FLIR92 

a. Advantages of FLIR92 

Since NVThermlP uses the same basic MRTD prediction theory with 
modihcations to improve on FLIR92 (NVESD 2001), familiarity would be the primary 
benefit of continuing with FLIR92. These modifications are discussed in the next 
sections, regarding the advantages of NVThermlP. If FLIR92 operates at a sufficiently 
accurate and precise level to satisfy the increasingly more stringent requirements in 
operational use of TAWS, then there is no reason to change. Also, NVThermlP 
does not calculate a MDTD, which is more useful for indicating thermal imager 
performance for point sources and aperiodic targets. TAWS uses the MDTD for 
certain calculations, when available. For these two reasons, if a simple modihcation 
to FLIR92 to include optical turbulence improves the FLIR92 MRTD prediction, then 
it would be advantageous to continue using FLIR92. 
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b. Disadvantages of FLIR92 

There are two major characteristics of FLIR92 that may lead to errors 
in performance predictions for staring arrays in sensors being modeled. Improved 
sensitivities in staring arrays mean that the limitations of the eye in detecting con¬ 
trasts is now restricting the performance of these sensors (NVESD 2001), but FLIR92 
MRTD calculations do not include any adjustments due to less-than-perfect contrast 
perceptions of the eye. Likewise, limitations on detector size, spacing, and hll factor 
(ratio of active cell area to total array area (Driggers 1999)) can cause under-sampled 
imagery for staring sensors. To avoid the under-sampled imagery, FLIR92 has an 
absolute cutoff at the half sample rates of imagers (Nyquist frequency), but when 
used with the Johnson criteria, it can lead to incorrectly negative predictions for 
most staring imagers (NVESD 2001). FLIR92 also does not account for atmospheric 
turbulence, which may be particularly signihcant in the highly turbulent regime just 
above a hot noon-time desert. 

4. Advantages and Disadvantages of NVThermlP 

a. Advantages of NVThermlP 

NVThermlP is considered an advance over FLIR92 primarily because it 
addresses the two possible errors in FLIR92 performance predictions regarding staring 
arrays, as discussed above. NVThermlP also does include atmospheric turbulence in 
a fairly rudimentary form. 

h. Disadvantages of NVThermlP 

There are three major hesitations in switching to use NVThermlP with 
TAWS. First is a lack of familiarity. Otherwise, since NVThermlP was developed to 
address weaknesses in FLIR92, it should be advantageous to use. For instance, NV¬ 
ThermlP uses the Targeted Task Performance (TTP) Metric instead of the Johnson 
Metric to provide a better performance estimate. This result is then used, like with 
FLIR92, to predict the target acquisition range performance for given sensors. Sec¬ 
ond, while NVThermlP predicts the MRTD, it is a laboratory MRTD that is not 
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comparable to the MRTD from FLIR92. Instead, NVThermlP predicts the CTF 
(NVESD 2001), which must be manipulated before insertion into TAWS. Finally, 
NVTHermlP does not produce an MDTD, but TAWS uses the MDTD for certain 
applications. 

5. Phenomena Missing from FLIR92 and NVThermlP 

Despite improvements in NVThermlP, there are details missing from both 
FLIR92 and NVThermlP. One important lack is the optical turbulence. NVThermlP 
does include an average optical transmission input, but neither model considers the 
effect of optical turbulence variations along the optical path from the target to the 
sensor. The next section discusses this in more detail. 

Neither FLIR92 nor NVThermlP have any adjustments for polarization of the 
target. Since the system response varies depending on the polarization, this lack will 
be especially apparent as imager sensitivities increase. Both models work best for 
nearly symmetric targets, but may poorly represent more realistic cases. 

B. MOTIVATION TO TEST IMPORTANCE OF OPTI¬ 
CAL TURBULENCE 

Recently, sensors with higher resolution have been developed. These IR im¬ 
agers are more sensitive to smaller targets, or a higher spatial frequency. In order 
to effectively model the responses of these sensors, it is becoming necessary to con¬ 
sider and perhaps include previously negligible effects like optical turbulence. At the 
sensitivities of the new sensors, the range for minimum detectability may be signif¬ 
icantly affected by slight variations in the optical turbulence along the path from 
the target to the sensor, perhaps even as regards a vertical pathlength through the 
atmosphere. Optical turbulence causes blurring around the edges of images (Fante 
1980) and makes it more difficult to resolve the image of the target well enough to 
identify it. According to the IREO Handbook Vol.3 (1993), the effects of pressure 
variations due to atmospheric turbulence are typically negligible compared to those 
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due to temperature and humidity fluctuations. In either very hot atmospheres such as 
those that may be operationally important in South Asia, or in very sensitive appli¬ 
cations - such as detection of a very small target, the effect of atmospheric turbulence 
may become significant. Recent improvements in sensor responsivities have increased 
the nnmber of operations in the sensitive applications category and therefore exam¬ 
ining the potential impact of optical turbnlence along the optical path throngh the 
atmosphere is crucial. 

C. PROJECT OVERVIEW 

Two comparisons to test possible improvements to the MRTD calculation were 
accomplished. 

First, a MatLab version of the operationally used FLIR92 was modified to in¬ 
clude variations in optical turbulence along the pathlength. Then the original FLIR92 
was compared to the modified FLIR92 to investigate the significance of improvements 
in the MRTD. Resultant MRTDs were inserted into TAWS and the maximum detec¬ 
tion ranges were compared. The generic IR sensor in TAWS was used and marine 
and desert environments were tested. 

Second, the MRTDs calculated by FLIR92 and by NVThermlP were com¬ 
pared. All input parameters were held constant for each model. The resnltant MRTDs 
from FLIR92 and NVThermlP were entered into TAWS. The atmospheric conditions 
were varied to test the marine and desert environments. In TAWS, the generic IR 
sensor was used for all tests, modified to include the MRTD from either FLIR92 or 
NVThermlP. 

1. Comparison Parameters 

To compare FLIR92 and NVThermlP nnder the conditions described, param¬ 
eters calculated by each were compared. Factors to compare and evalnate in order to 
determine the significance of any differences are the horizontal and vertical MRTDs 
and system MTFs at specified spatial frequencies. Again, the extrema of frequencies 
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are the regions of particular interest in this study. Because NVThermlP outputs the 
CTF, in order to compare, the resultant CTF was multiplied by 2 times the scene 
contrast temperature (SCT) to obtain the MRTD (Jon Hixson, personal communica¬ 
tion). 

The terms ‘horizontal’ and ‘vertical’ directions as used in this report refer to 
the along-bar and cross-bar directions of a four-bar target, as shown in Figure 1. 



h axis 

Horizontal Target 
Orientation 



Vertical Target 
Orientation 


Figure 1. Horizontal and Vertical Directions with Respect to a 4-Bar Target. 


2. Comparisons in TAWS 

Once the MRTDs were determined from FLIR92 and NVThermlP, the re¬ 
sults were entered into TAWS. Then, the range output from TAWS was compared 
for the various different iterations to determine how much and how signihcant the 
improvements to the maximum detection range in TAWS were. 
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II. THERMAL IMAGING SYSTEMS 
PERFORMANCE MODELS 

A. FLIR92 

1. MTFs Included in FLIR92 

FLIR92 uses a total system MTF to calculate the MRTD. The total system 
MTF is the product of the prehlter and posthlter MTFs, including the noise hlter 
MTFs. For most component MTFs, FLIR92 can supply a model prediction, but 
allows for a user to directly input a measured or previously known MTF for the 
component. Generally, the measured MTFs improve performance estimates and are 
preferred over the MTFs predicted and calculated by FLIR92. 

In FLIR92, the horizontal and vertical prehlter MTFs for a sensor are the 
products of the optical, detector spatial, focal plane array (FPA) integration, and 
image motion MTFs in the horizontal and vertical directions respectively. FLIR92 
allows user-dehned MTFs for all but the FPA integration MTF. 

For scanned observing systems, hltering in the temporal domain must be trans¬ 
lated into the spatial domain. The detector temporal MTF, electronics high and low 
frequency response MTFs, and boosting MTF are included in the temporal posthlter 
MTF that describes this translation. The temporal posthlter MTF in FLIR92 only 
applies for scanning systems and is only dehned in the horizontal direction. FLIR92 
allows a user-dehned temporal posthlter MTF. 

Both scanning and staring systems have spatial posthlters that model the 
system processing components. The spatial posthlter MTFs are the electro-optical 
multiplexer MTF, the digital hlter MTF, the display MTF, the Charge-Coupled De¬ 
vice (CCD) charge transfer efficiency MTF, the display sample and hold MTF, and 
the eye MTF. All are dehned in both the vertical and horizontal direction. 

FLIR92 allows for an extra user-determined prehlter MTF and three extra 
posthlter MTFs to enable adequate description of the system. All the following 
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MTFs are standard in FLIR92 and are summarized from the Analyst’s Reference 
Guide (1993) and compared against original sources or the IREO Handbook (1993) 
as specihed. For scanning systems, the scan direction is assumed to be horizontal to 
the bars of the target. 

a. Optical MTF 

The optical MTF includes the diffraction-limited MTF and geometric 
blur MTF. When available, the optical MTF may be replaced by a user-dehned 
version determined from direct measurements or ray tracing program predictions. 
The diffraction-limited MTF describes the resolution limitations due to diffraction in 
the optics of the sensor and is valid for a system with a circular, clear aperture {IREO 
Handbook Vol.4 1993). It it of the form 


HMs) 


2 

TT 


/ 

arccos 

V 


(Ml 

[Do 



( 2 . 1 ) 


where A is the wavelength for diffraction in jam and Dq is the optics aperture diameter 
in mm. Since arccos(a;) is only physical when —1 < a; < 1, then the spatial frequency 
fs, in cycles/mrad, must satisfy fs < The diffraction limited MTF in the hori¬ 
zontal direction is illustrated in Figure 2. For this report, the same inputs are used 
in the horizontal and vertical directions, so the MTF in the vertical direction is the 
same. 

Apertures that are partially obscured or are not perfectly circular, or 
aberrations in the optics may cause blurring of the image (Goodman 1968). The 
optics MTF includes a geometric blur term to describe this phenomenon. If there is 
no user input for the geometric blur MTF, then FLIR92 assumes the blur is Gaussian 
and the MTF has the form 


ffgb(fs) = exp(-27rVV.^)- 


(2.2a) 
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a, in mrad, is the standard deviation for a circular Gaussian blur distribution. The 
IREO Handbook Vol.4 (1993) clarihes: is 


a 


2 



(2.2b) 


where w is the Gaussian blur spot diameter in mrad at the ^ point. 

Figure 3 shows the geometric blur MTF in the horizontal direction. 
Again, with the inputs used in this report, the MTF is identical in the vertical direc¬ 
tion. 

b. Detector Spatial MTF 

Also included in the prehlter MTF is the detector spatial MTF that 
compensates for the hnite size of the detector. FL1R92 assumes a rectangular detec¬ 
tor geometry. Like for the optical MTF, the IREO Handbook Vol.5 emphasizes the 
importance of a user-specihed detector spatial MTF, particularly if the rectangular 
detector geometry approximation is grossly inaccurate. 


Diffraction MTF H 



Figure 2. Diffraction Limited MTF in the Horizontal Direction. 
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A Signal Processing in the Element (SPRITE) system, incorporating 
the signal time delay and integration into one elongated detector element {IREO 
Handbook Voh3 1993) has a horizontal detector spatial MTF of the form 

sm{7ilhfs) f I \ 


Hdsifs) = 


T^hfs \l + {2'KQfsY^ 


(2.3a) 


where Q is the diffnsion length and Ih is the instantaneons held of view (IFOV) of 
the readont length, both in mrad. The IFOV is the solid angle throngh which a 
detector is sensitive to radiation and in a scanning system. This is the snbtended 
solid angle when all scanning motion stops. The vertical detector spatial MTF for 
SPRITE systems and both horizontal and vertical MTFs for non-SPRITE scanning 
and staring systems is 


Hds{fs) = 


sin(7rh^/^ 


(2.3b) 


where 5z is the detector IFOV in mrad (Boreman 1990). The IREO Handbook Vol.4 
(1993) also points ont that this form is only valid when the point spread fnnction of 
the detector is a rectangnlar fnnction in direct space, bnt it is a reasonable approxi- 


Blur MTF H 



Fignre 3. Geometric Blur MTF in the Horizontal Direction. 
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mation in most cases. The MatLab version of FLIR92 used in this analysis assumes 
a non-SPRITE system and uses Equation 2.3b. Figure 4 shows the form of the de¬ 
tector spatial MTF in the horizontal direction, and again it is identical in the vertical 
direction. 


Det aprtr MTF H 



Figure 4. Detector Spatial MTF in the Horizontal Direction. 
This plot is assuming a non-SPRITE scanning or staring system. 


c. Focal Plane Array Integration MTF 

For scanning systems, the prefilter MTF also includes a focal plane 
array integration MTF introduced in the horizontal direction by the finite integration 
time of the detector. If is the detector integration time in s and Vg is the scan 
velocity in mrad/s, the MTF has the form of diffraction, or 


HMs) 


sm^TifsVgti) 


(2.4) 


This MTF is only signihcant for scanning systems because staring systems do not 
require compensation for hnite integration times since they do not move during de¬ 
tection. 
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d. Sample-scene Phase MTF 

When considering sampled systems, the target location on the sampling 
grid can cause a phase effect, or aliasing effect. The aliasing of high frequency noise 
into the bandpass of the sampled signal can produce false signals or greatly increased 
noise or interference {IREO Handbook Voh3 1993). The MTF representing this phase 
effect is 


Hssp^fs 



(2.5a) 


Here, tat is the Nyquist frequency in cycles/mrad and is the phase angle (rad) be¬ 
tween the target under consideration and the detectors operating at rj^. The Nyquist 
frequency, or Nyquist limit, is the frequency at which no useful information is trans¬ 
mitted and is taken as half the scene sample frequency. It is given as 



(2.5b) 


where s^ is the number of samples per detector IFOV. Input frequencies above the 
Nyquist limit are likely to appear as aliasing signals at the lower frequencies {IREO 
Handbook). At optimal operation, ©^ is set equal to 0, but for average conditions ©^ 
is set equal to 0.785 rad (45°). In the inputs used for the MatLab version of FLIR92, 
©^ is set to 0, or for optimal operation. This results in an MTF in the horizontal and 
the vertical direction that is constant at unity for all frequencies since cos(O) is unity. 

e. Image Motion MTF 

Image motion MTFs are included in FLIR92 to account for how the 
thermal system moves with respect to the scene being imaged. Two causes of these 
movements may be linear motion of the sensor system or vibration of the sensor 
platform. These movements are represented in the linear image motion MTF, random 
image motion MTF, and sinusoidal image motion MTF. 

The linear motion MTF explains smearing of the scene across the detec¬ 
tors due to movement of the system platform. How signihcant the smear is depends 
on how fast the platform is moving and how long the detectors are exposed to the 
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scene during platform movement. Given Vr is the relative image velocity in mrad/s, 
the linear image motion MTF may be described by 


Hmlifs) = 


( 2 . 6 ) 


TTVrtifs 

When ti is not specihed, the integration time defaults to the held time. Typically, the 
linear image motion is only signihcant for staring systems because scanning systems 
expect the sensor to be moving, or scanning, during detection. The MatLab version 
of FLIR92 sets this MTF to unity for all frequencies. 

For simplicity, the random image motion is assumed to be Gaussian. 
The IREO Handbook Vol.4 (1993) notes that this Gaussian term includes the random 
vibration of the sensor with respect to the scene, uncertainty in the sensor position, 
and other causes of image jitter. The associated MTF is 

HmrU's) = exp{-27r^ alosfs) (2-7) 


where <j\os standard deviation of the total random vibration from all sources, 

in mrad. Figure 5 illustrates the form of the random motion MTF in the horizontal 
direction. The vertical direction is identical with the given inputs. 

Finally, the sinusoidal vibration of the platform may be described by a 
zero-order Bessel function, as shown in the IREO Handbook Vol.5 (1993), 




( 2 . 8 ) 


A is the amplitude of the assumed sinusoidal vibration in mrad. The zero-order Bessel 
function has the form shown in Fig 6, similar to a damped cosine function. Using the 
given inputs, though, A is set to zero and H^s is thus held constant at unity for all 

frequencies since Jo(0) = 1. 

/. Detector Temporal MTF 

Because the detector has a hnite response time, the system has a low- 
pass temporal hlter given by 

(2,9) 


HmUi) = 
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where ft is the temporal frequency in Hz and is related to the spatial frequency by 
the scan velocity since ft = fsVs- fdt is the 3dB frequency in Hz of the detector 
temporal response. The IREO Handbook Vol.4 further elaborates on Equation 2.9 - 
the hnite temporal response to radiation incident on the detectors causes blurring in 
the scanning direction, which is why H^t is included in the horizontal for scanning 
systems. 

g. Electronics Frequency Responses 

The frequency responses for systems that are AC-coupled are given by 
multiple pole RC hlters. The IREO Handbook (1993) gives the single-RC low-pass 
and high-pass hlter MTFs, which are then raised to the nth power, where n is the 
number of hlter poles, to arrive at Equations 2.10a and 2.10b. 

For the low frequency response, the high pass hlter MTF is 


Hehpift) 



(2.10a) 


Motion Rnd MTF H 



Figure 5. Random Motion MTF in the Horizontal Direction. 
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where fehp is the electronics 3dB frequency in Hz. The low frequency response is only 
signihcant at extremely low frequencies, so FLIR92 does not include it in the total 
MTF, but Hehp is used to calculate noise bandwidths. 


For the high frequency response, the low pass hlter MTF is 



(2.10b) 


where like in the high pass hlter case, feip is the 3dB frequency in Hz and n is the 
number of hlter poles. This high frequency response MTF is included in the total 
MTF since the high frequency response is signihcant at a greater range of frequencies. 

h. Boosting MTF 

For the scanning system, an aperture correction MTF, or electronic 
boosting MTF, is required. The boost emphasizes higher frequencies to compensate 
for the reduced depth of modulation that typically occurs at higher frequencies due to 
the less-than-ideal aperture response {IREO Handbook Vol.4 1993). It has the form 



( 2 . 11 ) 


Zeroth Order Bessel Function 
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Figure 6. Zeroth Order Bessel Function 
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fb is the frequency in Hz of the maximum boost and Ba is the amplitude of the boost 
at a maximum frequency fmax- As the IREO Handbook Vol.4 observes, the boost 

MTF is applied along the scan direction. 

i. Electro-Optical Multiplexer MTF 

For the LED electro-optical multiplexer, the MTF is of the form of 


diffraction. 


kdeomi^fs^ 


sm{nfsSLED) 


( 2 . 12 ) 

eom is not included 


'^fs^LED 

and 6led is the angular subtense of the LED element in mrad. H, 
in the MatLab version of FLIR92. 

j. Digital Filter MTF 

Digital Liters have a linear phase symmetrical impulse response which 
varies depending on whether there is an even or an odd number of samples (iV). Eor 
N odd, 

2nif, 


(N-l)/2 

kddfi^fs^odd 'y ' Oj COS 

i=0 


fa 


For N even. 




N/2 

= COS 

i=l 


27r {i - I) fs 


fa 


(2.13a) 


(2.13b) 


For these MTFs, Oj is the digital filter coefficient and fco is the filter cnt off frequency 
in cycles/mrad. In the MatLab FLIR92, this MTF is set to unity for all frequencies 
in both the horizontal and vertical directions. 

k. Display MTF 

Althongh many recent thermal imagers use a cathode ray tube (CRT) 
display to communicate collected information to the observer, if a non-CRT display is 
being used, FLIR92 requires a user-specified MTF. Otherwise, the model calculates 
the total MTF based on a CRT display. In FLIR92, the phosphor spot Inminance 
intensity is assnmed to have a Gaussian distribution with a a valne either specified 
by the user or calculated by the model from the relationship 


a = 


\ 


log(0.025) 
■ 27r2/i 


(2.14) 
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with the raster frequency, /at^ in lines/mrad, 



(2.15) 


Here the model uses the TV limiting resolution criterion that the MTF equal 0.025 at 
the raster frequency. Nr is the number of lines on the CRT that are active. (3 is the 
vertical FOV in mrad. Using these quantities, the display MTF may be expressed as 



(2.16) 


Since this result is only valid for the specihed resolution criteria, modihcations for 
other resolution criteria may be made by multiplying a calculated a by an appropriate 
conversion factor and then entering the result directly into the model. The IREO 
Handbook Vol.5 (1993) again points out that a follows Equation 2.2b. In the FLIR92 
MatLab version used in this report, a was hard coded to 0.01 mrad, with the resultant 
MTF shown in Figure 7. 
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Figure 7. Display MTF. 

The standard deviation of the Gaussian distribution is set to 0.01 mrad. 
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1. CCD Charge Transfer Efficiency MTF 

For a CCD system, Sequin and Thompsett (1975) explained the charge 
transfer efficiency MTF as due to the non-unity charge transfer efficiency of the 
system. The CCD charge transfer efficiency MTF depends on the nnmber of gates 
in the transfer from the detector to the outpnt ampliher [N) and the charge transfer 
efficiency at each gate (e) along with the Nyqnist frequency, or the sampling freqnency 
of the strnctnre. As indicated in the IREO Handbook Vol.4, it may be expressed as 

Hccoifs) = exp(^-N{l - e) 1 - cos j (2.17) 


where r^r is again the Nyqnist freqnency in cycles/mrad (see Section II.A.l.d). N 
may be calculated in the model by assnming that the system has a interline transfer 
scheme, bnt often is specihed by the user. In the MatLab program being used, the 
CCD Charge Transfer Efficiency MTF was set to nnity for all freqnencies, in both 
the vertical and horizontal directions. 

m. Display Sample and Hold MTF 

Given that Sg is the sampling apertnre in mrad, the display sample and 


hold MTF has the form 


HdshUs) = 


sm{7rSsfs 


(2.18) 


Fignre 8 shows the horizontal display sample and hold MTF, bnt the vertical is set 
in the MatLab code to be nnity for all freqnencies. 

n. Eye MTF 

The eye MTF encapsulates the influences of the observer. In FLIR92, 
the eye MTF is based on work by Kornfeld and Lawson (1971). The observer may be 
able to improve and optimize the system display by adjusting the gain and level, the 
viewing distance, and the like. In this case, the eye MTF is considered non-limiting 
and is considered a constant of nnity since no degradation of the spatial freqnency 
response would be expected. In other words. 


HeyeUs) = 1 - 0 . 


(2.19a) 
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On the other hand, the observer may not be able to improve or adjust the display 
to peak performance. If this is the case, the eye limits the clarity of the image, so 
there is a degradation in the spatial frequency response. The MTF to describe this 
degradation, the limiting eye MTF, is of exponential form, 

Heyeifs) = exp • (2.19b) 

Here, M is the system magnihcation and F is the spread function width of the eye in 
mrad. F depends on the ambient light level (Kornfeld and Lawson 1971). Since the 
observer is often preoccupied or unable to adjust the display, FLIR92 defaults to the 
limiting eye MTF. In the MatLab code, F is hard coded to 2.196 mrad. The resultant 
MTF is shown in Figure 9. 

o. System Noise Filter MTF 

The system noise hlter MTF includes only the detector noise and is 
composed of all the noise hlter MTFs. These MTFs - Hdt, H^hp, roll up of the temporal 
posthlter MTF (Htpf), and the spatial posthlter roll up {Hspf) are discussed above. 


Sample hold MTF H 



Figure 8. Display Sample and Hold MTF in the Horizontal Direction. 
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The system noise filter MTF is not included in the total system MTF (Hsys) since 
the component MTFs are already included in the total MTF, but will be referenced 
in future sections. It may be expressed as 

HMfs) = Hdi{fs)Hehp{fs)HTPF{fs)HsPF{fs) (2.20) 

where for staring systems, HTPFifs) is unity. The horizontal system noise MTF is 
shown in Figure 10. 

p. Total System MTF 

The total system MTF is the product of all the component MTFs de¬ 
scribed above. For the generic case described in this report, the total horizontal MTF 
is shown in Figure 11. The downward curve of the MTF illustrates how much the 
image is degraded, particularly at higher spatial frequencies, corresponding to small 
targets. 


Display brt eye response 



Figure 9. Eye/Brain Response MTF. 
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2. Sampling 

a. Limits to Defined MRTD 

In FLIR92, the MRTD is only defined for a periodic target such as a 
4-bar target. The target must have a 7:1 aspect ratio, as shown in Figure 1 and the 
four bars should be fully resolved by the observer in order to specify an MRTD. In 
thermal imagers, the cutoff frequency for the observer to be able to fully resolve the 
target s the system’s Nyquist frequency. By definition, then, no MRTD can be given. 
This limits FLIR92 to MRTD prediction below the Nyquist limit. 

Since many targets are aperiodic, this artificial method of limiting the 
MRTD prediction may be pessimistic. Information relating to the MRTD may be 
available above the Nyquist limit, but methods for obtaining and quantifying this 
information are not available for use in FLIR92. 


Horizontal Noise MTF 



Figure 10. Horizontal System Noise MTF. 
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b. Sampling Rate and Lattice Spacing 

The IREO Handbook Vol.4 (1993) explains that imagers sample spa¬ 
tially and temporally. Spatially, the thermal imagers sample at a rate of 

R. = ( 2 . 21 ) 

Oz 

This is half the Nyqnist freqnency, as can be seen by comparing Eqnations 2.5b and 
2.21. The sample lattice spacing is then the inverse of Eqnation 2.21, FLIR92 
actnally calcnlates the effective sampling rate and effective lattice spacing by mnlti- 
plying together the horizontal and vertical components {z = h and z = v) and then 
taking the sqnare root. 

c. Optimizing Scene Phasing 

As discnssed in Section II.A.l.d, the MRTD measnrement may be de¬ 
graded by phasing between the target and the detector grid. There is an MTF 
inclnded to acconnt for this phase difference, particnlarly when dealing with staring 


Horizontal System MTF 



Fignre 11. Horizontal System MTF. 
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systems where the MRTD degradation is pronounced when there is a misalignment 
between respective targets and detectors under consideration. 

3. System Noise 

FLIR92 characterizes second-generation thermal imaging systems and as such, 
the system noise in MRTD predictions is modeled using a scaling factor that multiplies 
the random spatio-temporal noise by the amount of excess system noise (Kennedy 
1990). To simplify the model, the system noise is reduced to components that add in 
quadrature. Appropriate eye spatial and temporal integration effects are considered 
for each noise component included. 

The Night Vision and Electronic Sensors Directorate (NVESD) of the US Army 
in 1990 established a theoretical framework and a standard laboratory measurement 
procedure to better characterize the noise in the system (Webb et.al. 1991). This 
method of noise analysis isolates the system noise into eight components. These 
components are listed in Table I along with possible sources of the noise. 


Table I. 3-D Noise Component Descriptions in FLIR92 
Subscripts of noise components shown in table indicate dimension: t temporal, v 
vertical spatial, h horizontal. Table derived from Analyst’s Reference Guide (1993). 


Noise 

Description 

Potential Source 

(^tvh 

random spatio-temporal noise 

basic detector temporal noise 

^tv 

temporal row noise 

line processing, j, read-out 

O'th 

temporal column noise 

scan effects 

^vh 

random spatial noise 

pixel processing, detector-to-detector 
non-uniformity 

(Jy 

hxed row noise 

detector-to-detector non-uniformity, | 

CTh 

hxed column noise 

scan effects, detector-to-detector non¬ 
uniformity 


frame-to-frame noise 

frame processing 

S 

mean of all components 



Here, atvh is the basic detector noise which is often characterized by the NETD 
(see Equation 2.22, next section). The NETD is a sensitivity parameter. It is dehned 
as the temperature difference between the target and the background that will produce 
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a peak signal to rms noise ratio of unity at the output of a reference filter. This is how 
the NETD gives a rough estimate of how noisy a sensor is compared to the signals 
being detected. All other components besides atvh must be added in quadrature 
in order to completely characterize the thermal imaging systems noise (Scott et. 
ah 1993). In the current operational FLIR92, atvh is predicted but the remaining 
components are measured or estimated. FLIR92 includes estimates of 3-D noise for 
generic scanning and staring sensors so that when measurements are not available, 
a calculation may still be made. More detailed descriptions of the system noise are 
given in the sections following. 

a. Defaults for 3-D Noise Components 

In a given system where 3-D noise measurements have not yet been 
made or are not available, FLIR92 provides a set of default values. These default 
values are only given for what FLIR92 considers significant noise components and 
depend on the predicted random spatial-temporal noise, atvh- As discussed in Section 
II.A.3, the system noise sources add in quadrature, so only the most signihcant noise 
components are defaulted to non-zero components since the others would scarcely 
affect the outcome of the noise calculation. These default values are from a database 
of system noise measurements that were carried out by NVESD starting in April 1990. 
The values for each system in the database are normalized to atvh and averaged with 
other systems in the same class in order to determine which are the dominant noise 
sources in that class and which are so small they may be neglected. 

In scanning systems, the signihcant noise components are the temporal 
row and the hxed row noises, atv and Scanning systems typically show a wide 
range of variation in noise levels, so three default values are provided for both atv and 
av- Table II shows the model’s default values for scanning systems. 

In staring systems, the signihcant noise component is the random spa¬ 
tial noise, avh- NVESD found a single default value for avh to be sufficiently repre¬ 
sentative. This default value is shown in Table III. 
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Table II. 3-D Noise Components For Scanning Systems 
Subscripts of noise components shown in table indicate dimension: t temporal, v 
vertical spatial, h horizontal. Table derived from the Analyst’s Reference Guide 
_(1993)._ 


Noise Term 

Low Noise Default 

Moderate Noise Default 

High Noise Default 

^tv 

0.25atvh 

O.lbatvh 

\.t)atvh 

(Jd 

0.25atvh 

O.lbatvh 

\.t)atvh 


Table III. 3-D Noise Component For Staring Systems 
Subscripts of noise components shown in table indicate dimension: t temporal, v ver¬ 
tical spatial, h horizontal. Table derived from the Analyst’s Reference Guide (1993). 


Noise Term 

Noise Default 

^vh 

O.AOatvh 


b. Random Spatio-temporal Noise, atvh 
Since atvh is related to the actual bandwidth of the system, it may be 
measured directly at the output port prior to display. However, as outlined in the 
IREO Handbook Vol.4, it may otherwise be determined through the relationship to 
the NETD, 

where A/tv is the equivalent noise bandwidth for the NETD. In this case. A/p is 
the actual noise bandwidth which is associated with the system electronics before 
display. In order to determine the equivalent noise bandwidth, given that S{u) is 
the normalized detector noise power spectrum and Href{v) is the standard NETD 
reference hlter, the relationship used is 

OO 

A/jv = / dr. (2.23) 

0 

The actual noise bandwidth differs depending on whether a scanning or staring system 
is under consideration. When considering scanning systems, the noise bandwidth at 
the system output port is of similar form to the equivalent noise bandwidth, 

OO 

A/p = J S{r)[HTPF{iy)?dr. (2.24) 

0 


( 2 . 22 ) 
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For staring systems, the actual noise bandwidth is given by 




(2.25) 


and ti is the FPA integration time. 

For all classes of systems, staring or scanning, the general form of the 
random spatio-temporal noise is given in Equation 2.26, 


^tvh 


A2 


(2.26) 


7rrov^/D*(A,300)^(A)dA 

Ai 


Here, fno is the optical f-number. A/p is the system noise bandwidth as dehned 
in Equation 2.24 or 2.25, Tq is the optical transmittance, and is the detector 
area in cm^. The partial derivative is the thermal derivative of Planck’s Law in 
W/cm^/sr/pm. The detector noise-limited spectral detectivity is Zi)*(A,300) in cm- 
Hz 5/W and includes only noise components from the temporal noise sources since 
the spatial noise source contributions are included in the system noise correction 
functions described in the following section, II.A.3.C. Also, FLIR92 does not make 
any adjustments for changes to detector responsivity due to cold shielding, so any 
tweaking there must be off-line. 

c. Noise Correction Functions 

Since each 3-D noise component listed in Table I is assumed to be 
statistically independent, the total system noise, filtered, may be written as the root 
sum square of the noise components: 


+ alE,Ek.(f.) + ’’lEtE^.U.) + 

+ <ylEiM.)- (2.27) 


Here, Et, Eh^{fs) and Ev^{fs) represent the eye or brain temporal and spatial integra¬ 
tion associated with the noise components. The frame-to-frame term cr^Et has been 



dropped because the frame-to-frame noise at is almost always so small comparatively 
that it is negligible. The orientation of the MRTD target under consideration is in¬ 
dicated by the subscript 2 :. In this case, the temporal integration may be expressed 
in terms of Fr, the system frame rate in Hz, te, the eye integration time in sec, and 
at, the temporal sample correlation length. 

Et = (2.28) 

a rTe 

In FLIR92, at is assumed to be unity in MRTD calculations. The spatial integra¬ 
tions are simplihed in the FLIR92 model and may be expressed in terms of sample 
correlation factors ah and a^, horizontal and vertical sampling rates Rh and in 
samples/mrad, and spatial integration limits LhXfs) L,uXfs) mrad“^. These 
spatial integration limits are approximately the width and height of the MRTD bar 
target. For staring systems, the samples in each direction are assumed independent 
so that ah and q;„ are unity. For scanning systems, samples in the scan direction 
may not be assumed independent due to the motion of the scanner, so ah may be 
greater than unity. Similarly, samples taken perpendicular to the scan direction may 
be assnmed independent since the motion of the scanner is cross-directional, so a^ 
is nnity. Although FLIR92 uses this simplihed form, the exact expressions for the 
horizontal and vertical eye/brain spatial integration are given in Appendix B. 

In order to determine the noise terms for the horizontal and vertical 


MRTD calculations, it is assnmed that only noise components in the direction of the 
MRTD target degrade the MRTD. For the MDTD, since target orientation does not 
affect the calculation, the noise correction fnnction is independent of direction. Given 


the noise correction fnnctions in the horizontal. 



and in the vertical. 



(2.29) 


(2.30) 


29 



the horizontal and vertical MRTD noise terms may be expressed. In the horizontal 
direction, the noise term is given as 


^^hifs) = (ytvhkh{fs)jEtE^^{f,)EhSfs 


(2.31) 


and in the vertical, 


^v{fs) = (^tvhkv{fs)jEtE^^{fs)Eh^{fs). 


(2.32) 


For the MDTD, the noise correction function has the form 


■'(/s) — 1 + 


2 2 
^vh 771 — 1 I 
Z^Et + 


l^[E„(/,)]-' + ^[E*E,(/,)]- 


[E,if,)]-^ + ^[EtE,if,)y 


4. Calculating the Predicted MRTD 

The basic form for the MRTD calculated in FLIR92 is of the form 


(2.33) 


MRTD^if,) = 


^SNRTHCTtvhk^fs 

HsYsAfs) 


lE,EnAfs)E,Afs 


(2.34) 


where again indicates the orientation of the MRTD target in either the horizontal 
or the vertical direction. The denominator HsYsXfs) is the total system MTF, or 
the product of all component MTFs. kz{fs) is the 3-D noise correction function, as 
given in Equations 2.29 and 2.30. The threshold SNR required to resolve the MRTD 
target is SNRth and the eye/brain integrals are as described in Section II.A.c and 
Appendix B. The IREO Handbook Vol.5 further outlines this method of determining 
the predicted MRTD. More details about the included variables in this general form 
of the MRTD calculation are given below. The system MRTD for the given input 
parameters is shown in Figure 12. 

a. Psychophysical Constants SNRth and te 
In FLIR92, the threshold SNR to resolve the MRTD target and the eye 
integration time, SNRth and te, may be adjusted to tweak an MRTD prediction to 
better match a set of measurements. NVESD does have recommended settings, as 
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given in Table IV. Biberman (1973) showed using psychophysical data that SNRth 
is a function of the target spatial frequency, but the NVESD recommended value of 
2.5 is a representative average for optimal observing conditions. 


Table IV. NVESD Recommended Settings for Psychophysical Constants 
Settings in FLIR92 may be adjusted from NVESD recommendations. 


Psychophysical Constant 

NVESD Recommended Setting 

Units 

SNRth 

2.5 

- 

Te 

0.1 

s 


te depends on the background luminance. Luminances corresponding 
to a 0.1 s eye integration time show agreement with display luminances that NVESD 
measured in perception experiments in 1988. These experiments were conducted 
under conditions similar to those used in MRTD measurements such as a darkened 
room and optimal viewing. Observers set the display luminance to an average of 
0.15 milli-Lamberts. Higher ambient light conditions typically correspond to greater 


Test horizontal and vertical MRT 



Figure 12. System MRTD showing the Horizontal and Vertical Results. 
The 2-D system MRTD is calculated by taking the geometric average of the 
horizontal and vertical results, so it would lie between the two. 
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display luminances and a faster integration time may be appropriate in these cases, 
such as in the systems used in the held. 

b. 2-D MRTD Calculation 

FLIR92 can calculate a 2-D MRTD by taking the geometric mean of 
the horizontal and vertical MRTDs. In taking the geometric mean, each component 
is weighted equally and the mean is with respect to the spatial frequency axis. This 
causes the 2-D MRTD to asymptote at the mean value of the vertical and horizontal 
cut-off spatial frequencies. These cut-off spatial frequencies are determined by either 
the Nyquist limits or MTF roll off. 

5. Calculating the Predicted MDTD 

The basic form for the MDTD calculated in FLIR92 is 


MDTDU.) = 

where the MDTD is independent of target orientation, so kMoxifs) is independent 
of orientation and is given by Equation 2.33. Qh{fs), Qv{fs), and E^^fs) are dehned 
below and are equivalent for both scanning and staring systems. Along the horizontal 
direction, the eye/brain spatial integration differs for staring and scanning systems, 
so each is stated below. At is the target area in mrad^ and is related to the spatial 


frequency, fs by 


At = — 


(2.36) 


Equation 2.36 assumes an isometric target, but is a poor approximation for many 
operation targets that are not so uniform. All other variables are dehned as for the 
MRTD calculation. 

The Qhifs) and Qv{fs) integrals both include the total system MTF, Hsys- 
If represents the orientation, either horizontal or vertical, then the Qz{fs) integral 


Qzifs) — / [HsYSzi^)Y 


(2.37) 
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Ev{fs), the vertical eye/brain spatial integration, is defined the same way in both the 
staring and the scanning systems and is of the form 


Evifs) - — f [HsYsA^W[HNFA^)y 


sm 


(f) 


TTiy 

fs 


dv. 


(2.38) 


Ehifs) is defined differently for staring and scanning systems. The coefficient in front 
of the integral changes. For scanning systems it is of the form, 


E,{f,) = ^ j S{u)[HsysMY[H^f^{u)Y 


sm 


TTl/ 

fs 


du. 


(2.39) 


For staring systems, it is similar to the vertical direction given in Equation 2.38, 


A 

Ehifs) ^ ~ J [^SYShi^)Y[HNFhi^)y 


sm 


(f) 


TTU 

fs 


dv. 


(2.40) 


6. MRTD and MDTD Temperature Dependence 

Since background temperature affects the random spatial-temporal noise fac¬ 
tor, (Ttvh^ both the MRTD and MDTD calculations are influenced by the background 
temperature. In FLIR92, the default background temperature is 300 K, but the model 
can scale atvh to any background temperature between 240 K and 330 K. Generally, 
though, predictions made with a background temperature of 300 K are acceptable for 
most FLIR92 applications. In fact, atvh is the only 3-D noise component that FLIR92 
allows to be scaled to variable background temperatures. Temperature scaling only 
takes into account the background temperature; there is no scaling to account for the 
temperature dependence of transmittance in optical materials or noise from photon 
emissions by variably temperatured sensor components. For details of the tempera¬ 
ture scaling of atvh, see Appendix C. 


7. Johnson Criteria and FLIR92 

The Johnson criteria uses the resolving power of an imager under consideration 
to provide an estimate for how good the imager will be at target acquisition. The 
resolving power of the imager in cycles/mrad is multiplied by the target size in mrad to 
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get the number of cycles on target. This number of cycles on target is then compared 
to a table of the number of cycles on target needed to detect, recognize, identify, or 
otherwise perform target acquisition tasks. 

B. NVTHERMIP 

Prior to the introduction of NVThermlP, NVTherm was developed to rectify 
known errors in FLIR92, particularly as regards staring arrays. The FLIR92 model 
did not take into account the contrast limitations of the observer’s eye. Also, FLIR92 
limitations on detector size, spacing, and £11 factor caused under-sampled imagery. 
The new model attempted to compensate for the sampling artifacts by imposing a 
cutoff at half the sample rate of the imager, but when this result was used with the 
Johnson metric, resultant predictions were pessimistic. Although the basic MRTD 
theory used in FLIR92 went into NVTherm and NVThermlP largely unchanged, NV¬ 
ThermlP reworked the eye MTF contribution to improve the pessimistic results from 
FLIR92. NVThermlP predicts the CTF of thermal imaging systems, but can also 
predict a laboratory measurement of the the MRTD. NVThermlP also has completely 
redone the eye MTF. In this report, the CTF is used for comparison to the MRTD 
so the new eye MTF theory will be included in the comparison. 

1. MTFs Included in NVThermlP 

The MTFs included in NVThermlP are predominantly the same as those in¬ 
cluded in FLIR92. As in FLIR92, in most cases a user-defined MTF may be entered 
and is encouraged whenever possible, for more realistic predictions. All the follow¬ 
ing MTFs are derived from the NVThermlP User Manual (NVESD 2005) unless 
otherwise stated. 

a. Optical MTF 

The optical MTF still includes the diffraction-limited MTF and geo¬ 
metric blur MTF, both of the same forms as Equations 2.1 and 2.2a. NVThermlP 
assumes the optical MTF is the same in the horizontal and vertical directions. 
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b. Detector Spatial MTF 

NVThermlP calculates the detector spatial MTF using Equation 2.3b 
and does not distinguish between SPRITE and non-SPRITE systems. A user-specihed 
MTF would be required to clarify. Hds{fs) is not assumed to be symmetric in the 
horizontal and vertical directions. 

c. Focal Plane Array Integration MTF 

The focal plane array integration MTF is not included in NVThermlP, 
but could be entered by the user as one of the custom MTFs if so desired. 

d. Sample-scene Phase MTF 

The sample-scene phase MTF is not included in NVThermlP, but could 
be entered by the user as one of the custom MTFs if so desired. 

e. Image Motion MTF 

In NVThermlP, the only image motion MTF included is the random 
motion MTF that describes the vibration of the sensor platform. As in FLIR92 
Equation 2.7, Hmrifs) is assumed to be Gaussian. It is not assumed symmetrical in 
both directions. User-dehned MTFs for linear and sinusoidal motion could be entered, 
or a total motion MTF dehned. 

/. Detector Temporal MTF 

The detector temporal MTF is not included in NVThermlP, but could 
be entered by the user as one of the custom MTFs if so desired. 

g. Electronics Frequency Responses 

NVThermlP calculates a low pass hlter MTF for scanning sensors, only 
in the horizontal direction. The MTF is the same as Equation 2.10b. No high pass 
hlter MTF is outputted. 

h. Boosting MTF 

The boosting MTF is not included in NVThermlP, but could be entered 
by the user as one of the custom MTFs if so desired. 
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i. Electro-Optical Multiplexer MTF 

In NVThermlP, the electro-optical multiplexer MTF is calculated by 
Equation 2.12 unless the user specihes an MTF. It varies in the horizontal and vertical 
directions due to the possible differences in the LED size. 

j. Digital Filter MTF 

The digital hlter MTFs in NVThermlP are similar to those in FLIR92, 
Equations 2.13a and 2.13b, but have a factor to allow for dither. The dither factor, 
D, is set to 1 when the user chooses not to include dither. The modihed equations 
are given below. For N odd the MTF is 

{N-l)/2 


Hdfi^fs)odd ^ ' 


Oj COS 


i=0 


and for N even. 




N/2 

= ^ Oj COS 
i=l 


‘2nifs 
Dfco 

27r {i - i) /, 
Dfco 


(2.41a) 


(2.41b) 


k. Display MTF 

Though the display MTF in NVThermlP assumes a CRT display for 
presenting information to the observer, it includes a factor Ezoom that is not in 
FLIR92, Equation 2.16. Ezoom, Ezoom-, is related to the magnification in the system 
parameters and this magnification is 

Q ... arctan 


M = 


eye 


FOV., 


= 2 - 


\2Dvd) 


FOK, 


(2.42) 


where Dh is the display height in cm, Dvd is the display viewing distance in cm, 
and FOVy is the vertical sensor FOVvin mrad. When Ezoom is used, it is assumed 
that only part of the full FOV will be seen on the display screen. Ezoom examples 
are given in Table V. Values for M may vary from 1/6*^ to 200, but typically are 
between 0.5 and 20. 

Equation 2.16 is modihed to include the Ezoom factor by modifying a: 


a = 


S^FOVz 

FzoomD H 


(2.43) 
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Table V. Ezoom Examples 

Portions of total FOV area shown on display screen when different Ezoom values are 
used. Values obtained from the NVThermlP User’s Manual (2005). 


Ezoom 

M Faetor 

Vert FOV Seen 

Horiz FOV Seen 

Tot FOV Area Seen 

none 

1 

1 

1 

1 

single 

2 

1/2 

1/2 

1/4 

double 

4 

1/4 

1/4 

1/16 


where is the physical width of the Gaussian display spot in cm and FOVz is the 
FOV in either the horizontal {z = h) or the vertical {z = v) direction. 

Also, unlike FLIR92, NVThermlP includes MTFs for other possible 
display types besides CRT, but these are not discussed here since this report focuses 
on comparing the two models using the same inputs. 

l. CCD Charge Transfer Efficiency MTF 

The CCD Charge Transfer Efficiency MTF is not included in NVThe¬ 
rmlP, but could be entered by the user as one of the custom MTFs if so desired. 

m. Display Sample and Hold MTF 

The display sample and hold MTF of Equation 2.18 is again the same 
in NVThermlP, but is only calculated for scanning systems. Also, it is only applied 
when calculating the horizontal MTF. 

n. Eye MTF 

NVThermlP handles the eye MTF very differently from FLIR92. NV¬ 
ThermlP considers the human eye point spread function as a combination of three 
factors: the optics, the retina and the tremor. This leads to an MTF that may be 
expressed as the product of the component MTFs from the eye optics H^o, the retina 
Hret, and the tremor Hfrem'- 


H^yeifs) = HUfs)Hret{fs)HtreUfs)- (2.44) 

Overington (1976) outlines the theory of the human eye MTF and identihes the three 
component MTFs above. Based on Overington’s work, the general forms for the eye 
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MTF are found, but in In NVThermlP, these forms for the eye optical, retinal, and 
tremor MTF are all empirical (NVESD 2001). The eye optics MTF is 


= exp 


/ 

r 

io\ 

_ 

43.69^ 


V 

fo 

) 


(2.45) 


where Mgys is the imaging system magnihcation, and fo and io are dehned below. 


fo = exp{3M3-0.0216Dl\og,o{Dp)) 


= 0.7155 + 


0.277 


(2.46a) 

(2.46b) 


The variable Dp is the pupil diameter in mm and is valid if one eye is used. If two 
eyes are used, then NVThermlP reduces the pupil diameter by 0.5mm. Dp is dehned 
as 


Dp = —9.011 + 13.23 exp 


logio(/i 


(2.46c) 


21.082 
\ / 

where fi is the number of foot-Lamberts at the eye from the display and is /l = g^. 
Ld is the display luminance in milli-Lamberts. The retina MTF is dehned as 

I 1.21N 


HretU's) = exp -0.375 


fs 


M. 


sys 


(2.47) 


The tremor MTF, or the MTF of the eye due to tremor, is 


Htremifs) = exp -0.4441 


fs 


M. 


sys 


(2.48) 


Clearly, then, the MTFs depend on the pupil diameter, which in turn depends on the 
light level. 

Figure 13 compares the eye MTFs in FLIR92 and NVThermlP. Clearly, 
NVThermlP has a more optimistic estimate of the image degradation due to the 
observer’s eye, which reduces the known error in the FLIR92 modeling of the human 
eye. 
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o. System Noise Filter MTF 

The system noise filter MTF is merely the roll-up of all the system 
noise. In NVThermlP, there is a further factor for system noise that is requested as a 
user-input. NVThermlP considers that each system has a different hxed pattern noise 
associated with variations in detector gain and level offset. There are three options 
in the NVThermlP input: None, Noise Factor, and 3-D Noise. If None is selected, 
then the system is considered ideal with no variations in gain and level offset among 
the detectors. For Noise Factor, experimentally-determined independent horizontal 
and vertical factors are multiplied by the horizontal and vertical noise CTFs. The 
3-D Noise is discussed in a later section about system noise. 

p. Interpolation MTF 

NVThermlP includes an MTF for interpolation, or the process of in¬ 
creasing the image size by inserting ‘filler’ pixels between the original pixels. Interpo¬ 
lation may be done either vertically or horizontally, and if interpolation is selected for 


Compare Eye MTF from FLIR92, NVThermlP 



Figure 13. Comparing the Eye MTF in FLIR92 and NVThermlP in the Horizontal 
Direction. 
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both directions, only one is applied at a time. The results are then combined. NV- 
ThermlP offers three different methods for interpolation, Pixel Replication, Bilinear 
Interpolation, and Bicubic Interpolation. In this report, since FLIR92 does not offer 

interpolation, no interpolation is selected. 

q. Optical Turbulence MTF 

Unlike FLIR92, operational NVThermlP already includes an MTF for 
optical turbulence. This turbulence MTF assumes that is constant along the 
pathlength and so is an average MTF. is a parameter that describes the optical 
turbulence in the atmosphere and more thoroughly dehned in Chapter III. The 
formulation is very similar to the Goodman results for constant given in Equation 
A.98 in Appendix A. The equation is given below, 

Hatifs) = exp (^-57.4:afiClX-h ^ jj (2-49) 

where a is a constant dehned as |, /^ is here dehned in cycles/rad, is still in m“i, 
and A, ; 2 , and D are all in m. 


2. Sampling 

a. Sample Spacing 

For an imaging system, the sample spacing quantihes the limits due to 
sampling. For staring systems, the sample spacing can be calculated in the horizontal 
and vertical directions by 

ass. = (2.50a) 

where indicates the target orientation, either horizontal or vertical. is the number 
of detectors in the direction, ass. is in mrad. For a scanning system, the vertical 
sample spacing is calculated as per Equation 2.50a, but in the horizontal direction 
the sample spacing is diherent. There is no sample spacing in a continuously scanning 
system, so assh = 0. A scanning system that samples has a sample spacing requires 


a user-input of Nhifov-, or samples per HIFOV. Then the sample spacing is 




DASh 

Nhifov 


(2.50b) 
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DASh is the detector angular subtense (DAS) in the horizontal direction, in other 
words, the detector width divided by the focal length of the collecting optics. 

b. Sample Frequency 

As explained in Section II.A.2.b, the sampling frequency is just the 
inverse of the sample spacing. Equations 2.50a or 2.50b. Thus, the sampling frequency 

is 


fsamp 




(2.51) 


Since the horizontal and vertical sample spacing may differ, the horizontal and vertical 
fsamp may also differ. Again the half sample frequency is the Nyquist frequency. 


3. System Noise 

As mentioned in Section II.B.l.o, the 3-D system noise is handled in a similar 
manner to FLIR92. The 3-D system noise components are dehned in Table VI, 
essentially the same as Table I. 


Table VI. 3-D Noise Component Descriptions in NVThermlP 
Subscripts of noise components shown in table indicate dimension: t temporal, v ver¬ 
tical spatial, h horizontal. Table derived from NVThermlP User’s Manual (NVESD 
2005). 


Noise 

Description 

Potential Source 

^tvh 

random spatio-temporal noise 

Basic Detector Temporal Noise 

^tv 

temporal row noise, line bounce 

Line Processing, |, read-out 

O'th 

temporal column noise, column 
bounce 

Scan Effects 

(^vh 

random spatial noise, bi¬ 
directional hxed pattern noise 

Pixel Processing, Detector-to-Detector 
Non-Uniformity j 

(Jy 

hxed row noise, line-to-line non¬ 
uniformity 

Detector-to-Detector Non-Uniformity 

(Jh 

hxed column noise, column-to- 
column non-uniformity 

Scan Effects, Detector-to-Detector 
Non-Uniformity 

(Jt 

frame-to-frame noise, frame 
bounce 

frame processing 

s 

mean of all noise components 
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As in FLIR92, these eight noise parameters are derived using directional av¬ 
erages. atvh is predicted and all other components are estimated based on historical 
databases of measurements. In Table VI, the subscript that is missing indicates which 
directions were averaged so, for example, atv was averaged in the horizontal and 
in both the temporal and horizontal. Again similar to FLIR92, only certain noise pa¬ 
rameters are considered important in scanning and staring systems. Table VII shows 
the key parameters for scanning systems, where the noise term is normalized using 
the random spatio-temporal noise 

Table VII. Noise Values for Scanning Systems in NVThermlP 
Subscripts of noise components shown in table indicate dimension: t temporal, v ver¬ 
tical spatial, h horizontal. Table derived from NVThermlP User’s Manual (NVESD 
2005). 


Noise Term 

Low 

Moderate 

High 

(^vh/(^tvh 

0 

0 

0 

^v/ ^tvh 

0.25 

0.75 

1 

^h! (^tvh 

0 

0 

0 


Table VIII shows the key parameters for staring systems, where the noise term 
is again normalized with atvh- Table VIII compares to Table III. 


Table VIII. Noise Values for Staring Systems in NVThermlP 
Subscripts of noise components shown in table indicate dimension: t temporal, v ver¬ 
tical spatial, h horizontal. Table derived from NVThermlP User’s Manual (NVESD 
2005). 


Noise Term 

Low 

Moderate 

High 

(^vh/^tvh 

0.2 

0.5 

1 - 2 

/ ^tvh 

0.2 

0.5 

1 - 2 

(^tvh 

0.2 

0.5 

1 - 2 
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a. Random Spatio-temporal Noise, Otvh 
The random spatio-temporal noise is calculated similarly in NVThe- 
rmlP as in FLIR92, but includes a factor to account for the number of detectors. 
Again, atvh is calculated assuming a background temperature of 300K and the rela¬ 
tionship, 


(^tvh — 

ToVA^/£>*(A,300),||j(A)rfA 

Ai 


(2.52) 


where N is the number of detectors and all other variables are dehned as for Equation 
2.26. Note that Equation 2.52 differs from Equation 2.26 by a factor of tt in the 
denominator. 


4. Calculating the Predicted MRTD 

NVThermlP outputs three laboratory measurements of the MRTD: high gain, 
low gain, and user input. The basic equation used in all cases is 


MRTD^U,) 


2SCNTMpCTFeye{fs) a^Tl,QH4fs)QV4fs) 

(A^^rifs) - ^pc(/.))^L ^ SCNLp 


(2.53) 


Here, SCNtmp is the scene contrast temperature that generates the average display 
luminance in K, CTF^yeifs) is the naked eye CTF, Ai,ar is the area of the bar in the 
target in cm^, A^pc is the area of the space between bars in the target in cm^, Sl is 
the unitless normalized laboratory detector responsivity, Tdet is the noise variance in 
K-mrad-sT and QH^ and QVz are the horizontal and vertical noise bandwidth for the 
horizontal system CTF and the vertical system CTF. The laboratory conditions can 
be strictly controlled, which is important because in this calculation when Abar{fs) = 
Aspdfs), the MRTD will be undehned. 

The laboratory MRTD in NVThermlP is difficult to compare directly to the 
MRTD in FLIR92 for many reasons. Two of particular concern in this study were 
the high, low, and user gain of NVThermlP - there is no equivalent in FLIR92 - and 
the fact that NVThermlP calculates the MRTD for laboratory conditions which are 
highly controlled and not equivalent to the conditions assumed in FLIR92. Since 
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the laboratory MRTD outputted by NVThermlP is not directly comparable to the 
MRTD from FLIR92, a work-around was proposed in conversations with Jon Hixson 
(Army Research Labs, Fort Belvoir, VA). The CTF outputted from NVThermlP can 
be manipulated to obtain an MRTD comparable to FLIR92 by multiplying 2 x the 
SCT. NVThermlP produces a system CTF in the horizontal and vertical directions 
following Equations 2.54a and 2.54b: 


CTF„,JfJ 

CTFyJU 


\ 


lC'TF„,Jf,Jf + 


CTF„ If. 

^^noise \«/ ‘ 




SCNi 


tmp 


\ 


[CTF,, .(/.„)] 2 + 


CTFv (F 


SCN, 


trap 


(2.54a) 


(2.54b) 


where the eye CTFs and the noise CTFs in the horizontal and vertical directions are 
dehned as 


CTF„,JfJ 

CTFv,Jf..) 



( CTF,y, 

( Th ) N 

y 


V Smag) 

(2.56a) 

\Md^spMTFH^yXfs,) ^ 


( CTF,y,{ 

' fs. ) \ 


^a^Tl,QH,,rQKer{U ( 

^ SmagJ 1 

(2.56b) 

^Md^spMTFv^yM's.) j 

CTF ( "l 

^eye {SmagJ 
MdispMTFn^y, ifsh) 



(2.56c) 

MdispMT Fv^y^ (/s„) 



(2.56d) 


The Mdisp is the display glare and a is a proportionality constant of 169.6 Hz 2 . 
SCNfmp is solicited from the user and considered a constant throughout calculations. 


5. Detector Cooling in NVThermlP 

NVThermlP allows for the selection of an uncooled detector. Selecting this 
option requires a user input of a performance measurement for the uncooled array 
in terms of the measured detector noise, detector frame rate, f-number, and optics 
transmission. Then, NVThermlP will calculate a peak D* and integration time for 
the uncooled sensor. 
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6. Targeting Task Performance (TTP) Metric and NV- 
ThermlP 

Unlike FLIR92 and previous version of of NVTherm, NVThermlP uses the 
Targeting Task Performance (TTP) metric to predict the probability of successful 
task performance. Since this study does not utilize the range performance predictions 
from NVThermlP, no further details are provided here, but for more information, see 
the NVThermlP User’s Manual (NVESD 2005). 
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III. OPTICAL TURBULENCE MODEL 

CHOICE 

Optical systems operating in a turbulent atmosphere experience broadening of 
the point spread function, which is best described using a specihc atmospheric MTF. 
For a collimated beam passing through the atmosphere, turbulence distorts the shape 
of the wavefront and causes variations of intensity along the wavefront so that when 
the beam is brought back into focus by some optical system, the image formed has 
been altered by the atmospheric affects. Driggers (1999) notes that the atmospheric 
MTF may be expressed as two separate MTFs created by two different effects in the 
atmosphere: scattering due to aerosols and blurring due to turbulence. In this study, 
the MTF describing scattering due to aerosols is neglected. 

This study uses the mathematical model of the atmospheric turbulence MTF 
developed by Goodman (1985). Various mathematical determinations of the atmo¬ 
spheric parameter are considered for use in the atmospheric turbulence MTF 
calculation, in order to most appropriately represent through the whole atmo¬ 
sphere. 

A. OPTICAL TURBULENCE AND ATMOSPHERIC 0 % 

According to the IREO Handbook Vol.2 (1993), turbulence in the atmosphere 
creates random variations in the atmosphere’s index of refraction. These irregularities 
distort the wavefronts that pass through them and thereby cause image distortion 
or blur in imaging systems. Although there are different geometries for describing 
turbulent systems, the applicable one in dealing with airborne and ground-based 
sensors is spherical wave propagation. In this case, the propagating light comes from 
sources that are in or near the turbulence, as in the imaging of objects in the turbulent 
atmosphere. The sensors receiving the light are also in the turbulent atmosphere. 
Although this study focuses on the airborne sensor situation where a target on the 
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surface is imaged by an over-flying aircraft, it is important to note that it is the 
turbulence at the sensor that most strongly affects the image quality. 

Operationally, the best estimate of atmospheric turbulence would be deter¬ 
mined from an analysis or forecast, but when a local estimate is not available, the 
atmospheric turbulence may be predicted using models that estimate the atmospheric 
structnre parameter as described below. The atmospheric structnre parameter 
describes atmospheric variations in the index of refraction. Essentially, fluctuations 
in the index of refraction along the path between the target and the sensor causes 
some IR radiation to be randomly bent from the path, resulting in a blurred image 
at the sensor. The more tnrbulent the atmosphere, the greater the blnrring. is 
best dehned by Equation 3.57, 

C^ix) = [n{x) — n{x + r)]r“3 (3.57) 

where n{x) is the index of refraction at a point x in the atmosphere and n{x + r) is the 
index of refraction some distance r away from x. The over-bar indicates an average 
over the representative part of the atmosphere by either averaging n over time at one 
location or space at a very short time (Goroch 1980). 

Atmospheric variations in the index of refraction are described by The 
generic form of how varies in the atmosphere is given by Friehe (1977) as 

Cl = (^79 X (G| + 0.113Gtq + 3.2 x lO'^Gj) (3.58) 

where P is the pressnre in mb, T is the temperature in K, G|. is the temperatnre 
strncture fnnction parameter, Ctq the temperatnre-water vapor parameter and Cq 
the water vapor parameter, all in K^/m“i. The IREO Handbook Vol.2 (1993) states 
that the dry-air, or G^, term dominates in most applications since generally the G-f’Q 
and Cq terms comprise no more than 2% of the total G^. Since G|. is typically the 
dominant term in G^, in this stndy it will be considered the only contributor. All 
three structure parameters vary depending on location in the boundary layer (Fairall 
1982). Althongh it will not be addressed in this study, it should be noted that in the 
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microwave part of the spectrum, depends much more strongly on the humidity 
and Ctq and Cq are large contributors to the overall (Fante 1980). Also, in 
rare situations in the IR part of the spectrum, Ctq may be more signihcant {IREO 
Handbook Vol.2 1993). 

1. Turbulence Near the Ground 

There are two primary ways that the ground impacts the movement of air and 
causes turbulence. First, the free air stream flowing along the ground experiences 
friction due to surface roughnesses, which causes wind shear. Second, the ground 
may serve as a source or a sink for thermal energy of the air. Given sunny conditions 
during the day, the ground will be a source of heat because the sun warms the ground 
to a temperature higher than the air above it. This leads to thermal convection 
and instability. At night, the ground will act as a heat sink by radiative cooling, 
resulting in a ground temperature colder than the air above it. These conditions 
are considered stable. When the air and the ground are at the same temperature, 
atmospheric conditions are considered neutral. These extreme fluctuations in ground 
temperature are observed over land, but over the ocean the temperature difference 
between night and day is much less. The variation in temperature differences was the 
primary motivation for comparing FLIR92 and NVThermlP over desert and marine 
locations. 

Fairall et. ah (1982) propose that in the surface layer of the atmosphere, the 
structure function parameters have the forms of 

= T^Z-lfiO (3.59a) 

GJ = QlZ-lAfiO (3.59b) 

Ctq = ttqT^Q^^Z sy/Af{^) (3.59c) 

where T* and Q* are temperature and humidity scaling parameters in K and g/m^ 
respectively, f{^) is a dimensionless function, and ttq is the temperature-humidity 
correction parameter. Fairall et. ah (1982) gives an estimate value of 0.8 for ttq in 
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unstable conditions such as those being studied here. A is a constant, taken to be 
0.6. Z is the height above the surface in m. 

The dimensionless function is connected to the Richardson number, so it 
accounts for varying conditions depending on the stability. Wyngaard et. al. (1971) 
determined the empirical form. 


/(O 

= 4.9(1-7.0O“5 

e<o 

(3.60a) 

fiO 

= 4.9 

e = o 

(3.60b) 

fiO 

= 4.9(1 + 2.750 

e>0 

(3.60c) 


where < 0 corresponds to unstable conditions (as under consideration here), = 0 
is neutral and .^ > 0 is stable. Wyngaard originally split /(.^) into two segments of 
.^ < 0 and ^ > 0, but the neutral case has been added here to show when is a 
constant. ^ itself is dehned as the height scaled by the Monin-Obukhov length scale. 


^ KgZ (t* + 0 . 61 Ta^) 


(3.61) 


L is the Monin-Obukhov length scale and is dehned as the height over the ground 
where the mechanically produced turbulence from vertical shear balances the dissi¬ 
pative effect of negative buoyancy. In other words, the Richardson number equals 
unity. L may be expressed as 


L 


T 


Kg (t* -F 0 . 61 Ta^) 


(3.62) 


where k is the unitless von Karman constant that Fairall et. al. approximate to 0.35, 
g is the gravitational acceleration on earth taken as 9.8 m/s^, Q* is the humidity 
scaling parameter in g/m^, p is the density of air in kg/m^, and u* is the frictional 
velocity in m/s. In this case, the density of air is taken as 1.3 kg/m^. Ta is the 
temperature of the ambient air in the region of interest, in K. 

If is taken to primarily depend on as was assumed above, then Equations 
3.59b and 3.59c may be neglected and their terms in Equation 3.58 dropped. Davidson 
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et. al. (1978) indicate that above the ocean and within lOm of the ocean surface, 
0“^ varies as but above that lOm, varies as Z~i. This is better seen by 

solving Equation 3.59a for the conditions given in Equations 3.60, as shown below 
where again ■C = f > 


Cl 

= 4.9T|Z-§ [1 - 7^]"5 

e<o 

(3.63a) 

c| 

= 4.9T|Z“3 

e = o 

(3.63b) 

cl 

= 4.9T|Z-§ (1 + 2.750 

e>0 

(3.63c) 


Substituting Equation 3.61into 3.63 for the three cases of unstable, neutral, and stable 
respectively, yields 


C| 

3 * 

(3.64a) 

C| 

= 4.9T|Z"i 

(3.64b) 

C| 

~ 13.5T|Z3L-h 

(3.64c) 


In unstable conditions, |.^| 3> 1 (Equation 3.64a), so (1 — 7^) ^ can be approximated 
_ 2 

as (—7.^) 3. In the stable case, |.^| -C 1 (Equation 3.64c), so (1 + 2.75.^) is approx¬ 
imately 2.75^. Other authors indicate that for |,^| 1, may even be taken as 

independent of Z. 

The above discussion applies over the ocean, but Hall (1977) observed that 
in unstable conditions near land surface, still decreases with height by Z~3. The 

4 

Z~3 height dependence frequently continues beyond 10m over land, but is still only 
valid in the surface layer. 

2. Turbulence Above the Surface Layer 

In the so-called mixed layer above the surface layer, Fairall et. al. (1982) offer 
dehnitions for C^, Cq, and Ctq where in all cases falls off as a function of 
They are very similar to Equation 3.60, so in this study. Equation 3.60 is used beyond 
the lowest surface layers. 
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Above the mixed layer, one possible model of behavior is the Hufnagel 
model. Goodman (1985) offers a possible analytic approximation to how varies in 
the vertical and references the work by Hufnagel and Stanley (1964). This approxi¬ 
mation is 

y s 10 

^^(Z) = 2.7 X 10"^® 3 < f — j e-^ + e-^ (3.65) 

where Z is the location along the flight path in km and < > is the mean value 

of the squared wind speed. is again in m“i. This approximation reasonably rep¬ 
resents the average variation of in the middle to upper atmosphere but poorly 
represents the boundary layer. Variations of C\ in the boundary layer differ exten¬ 
sively depending on the stability of the boundary layer, as described above. Fante 
(1980) and other sources indicate that the Hufnagel approximation agrees fairly well 
with observation starting around 5 km above the surface. 

3. Turbulence Through Whole Atmosphere 

In this study, the atmospheric structure parameter G^ is taken to follow the 
bulk method below 200 m and the Hufnagel formulation above the 5 km point. It is 
worth noting that the Hufnagel formulation is an early numerical model of the upper 
atmospheric structure parameter and is generally only considered valid for above 5 
km in the atmosphere. For simplicity in MatLab programming, between 200 m along 
the pathlength and 5 km, an average value of 1 x 10“^® m“3 was chosen. It is near 
the surface that C\ varies the most and is the largest, so the actual choice for C\ 
above the surface layer was of less concern than that in the surface layer. The C\ 
calculation is summarized in Equation 3.66, 


G2 

= BulkMethod 

Z < 200m 

(3.66a) 

G2 

= lx 10“^®m“5 

200m < Z < 5000m 

(3.66b) 

G2 

= Hufnagel Method 

Z > 5000m. 

(3.66c) 


In order to test the validity of these assumptions, the bulk method and the Hufnagel 
formulation were plotted for several temperature differences. In Figure 14, the surface 
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temperature is taken to be 28°C and the ambient air temperature 25°C. Other param¬ 
eters are given in Table IX, where for this test case, generic parameters are chosen. 
The two methods intersect below 200 m, but at 200 m the values are on the order 
of 1 X 10 m 3, which agrees fairly well with results from Fairall et. ah (1982) and 
other sources. Values differed by less than an order of magnitude using the Hufnagel 
model at 5 km. The choice of a constant in the middle atmosphere simplihed the 
MatLab programming, but the small scale of and the relatively small difference 
between the two models around the 500 m to 5000 m layer suggested that choosing a 
constant was reasonably representative of the atmospheric conditions. Note that 
for simplicity in the MatLab programming, the bulk method, which describes C\ in 
the surface layer, is taken to be valid up to 200 m. This does not mean that the 
surface layer itself extends up to 200 m. Rather, the bulk method is assumed valid 
up to 200 m, so the bulk method may be assumed valid beyond the surface layer. 

Table IX. Parameters for C\ Test Cases. 

Generic parameters chosen to test the Bulk and Hufnagel formulations of C\. 


Parameter 

Value 

Units 

Rel Humidity Ambient Air 

0.8 

fraction 

Rel Humidity Sfc Air 

1.0 

fraction 

Wind Speed 

7 

m/s 

Sfc Pressure 

1012 

hPa 

Thermal Sfc Roughness 

0.02 

m 

Momentum Sfc Roughness 

0.02 

m 

Total Pathlength 

10,000 

m 


After comparing the bulk method and the Hufnagel method for a sample 
test case, the two methods were compared for conditions at White Sands and Point 
Conception (data give in Chapter IV). As indicated in Figures 15 and 16, the assumed 
constant value of 1 x 10“^® m“i appears fairly representative between the two models. 
Also, the exponential behavior of the bulk method near the surface, in both cases, 
indicates how much stronger turbulence is near the ground. 


53 




Note that while the surface bulk model and Hufnagel model were used this 
study, there are numerous other models of the atmospheric structure parameter. It 
is also important to note that this study does not consider possible surface inversions 
and neglects the stable case. 


B. MTF FOR OPTICAL TURBULENCE 

The Goodman (1985) formulation for atmospheric optical turbulence was cho¬ 
sen for this study. It can be expressed for the long-exposure and short-exposure case, 
although the focus here is on the long-exposure case. 

The long-exposure and short-exposure atmospheric turbulence MTFs differ 


Comparing Models Sfc T = 28 C, Amb T = 25°C 



10 


Figure 14. Comparison of the Bulk and Hufnagel Methods with Sfc Temp 28°C. 
At 28°C, the two methods intersect at 158 m. At 200 m, the bulk method has a 
value of approximately 1 x 10“^® m“3. The bulk method is used below 200 m. 
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substantially. In the long-exposure case, the exposure time is such that the phase and 
log-amplitude vary over the course of the exposure. Thus during a long exposure, the 
image recorded will be spread due to random variations in the tilt of the wavefront. 
To consider the short-exposure case, a random factor associated with the wavefront 
tilt is extracted from the MTF before the average is taken. This is not done for the 
long-exposure case. For the extremely short exposures, there will be no impact from 
the wavefront tilt and it is ignored in determining the MTF (Fried 1966). 

The derivation for the Goodman atmospheric MTF is summarized in Ap¬ 
pendix A. Using the generic atmospheric input parameters given in Table IX, the 
Goodman long exposure atmospheric turbulence MTF is orders of magnitude smaller 
than the other MTFs included in FLIR92. In Figure 17, Goodman MTF is subtracted 


Comparing Models Sfc T = 75 C, Amb T = 35 C 
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Figure 15. Gomparison of the Bulk and Hufnagel Methods at White Sands. 
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from unity to better illustrate the form of the MTF. The used in the Goodman 
calculations may be calculated using Equation 3.66. 

From the above discussions, it appears that is most signihcant near the 
surface. In one illustration of this, the Goodman atmospheric MTF is examined 
for a variety of slant angles starting with a path perpendicular to the surface, and 
ending with a virtually horizontal path through the surface layer, parallel to the 
surface. Figure 18 illustrates the Goodman MTF at various slant angles, for a constant 
pathlength of 100 m. Plotted is unity minus the Goodman MTF, as in Figure 17, to 
best observe this small scale MTF. Although Figure 18 makes it clear that even close 
to parallel to the surface, the Goodman MTF is small scale, still the resultant MTF 
is more signihcant at smaller slant angles. 


Comparing Models Sfc T = 31 °C, Amb T = 25°C 



Figure 16. Gomparison of the Bulk and Hufnagel Methods at Pt Gonception. 
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Figure 17. Goodman Long Exposure Atmospheric Turbulence MTF. 

Note the small scale of the atmospheric turbulence MTF. Plot is of (1 - Goodman 

MTF). 
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MTF 



Figure 18. Comparison of the Goodman MTF at Different Slant Angles. 
Perpendicular to the surface, the Goodman MTF is least significant and parallel to 
the surface, most signihcant. In all cases, the Goodman MTF is small scale. 









IV. 


METHODS AND DATA SETS 


A. SELECTION OF DATA LOCATIONS 

This study compares MRTDs for marine and desert environments, with and 
without an included MTF for atmospheric turbulence. Since this study is only an 
initial comparison, only two sites were selected for analysis. Because of access to ac¬ 
curate data, the White Sands Missile Range in New Mexico was selected for the desert 
environment and Buoy 46063, just off Point Conception near Santa Barbara, CA was 
selected for the marine environment. Althongh the two locations differ in nearly all 
atmospheric conditions, the critical difference is between the surface temperature and 
the ambient air temperatnre. As explained in Section III.A, the temperatnre struc¬ 
ture parameter is the largest contributor to variations in the atmospheric strnctnre 
parameter. 

B. SELECTION OF ATMOSPHERIC DATA 

Atmospheric data for inclusion in TAWS and for calculation of parameters 
in FLIR92 and NVThermlP were obtained from several sonrces. For the land lo¬ 
cations, atmospheric data was obtained from the Weather Underground website, 
www.wnnderground.com. The website provided all the needed TAWS inputs. An 
extreme example of the snrface temperatnre and ambient air temperature at White 
Sands was provided by Rick Sharkey and cohorts at the White Sands facility. These 
input parameters are given in Table X. 

For the atmospheric data over marine locations, observations from ships and 
bnoys were used. The team at AFWA provided ship data and npper air analysis for 
the Point Conception region. The inpnt parameters are given in Table XI 

In both cases, data was obtained for 27 Aug 07, at 2030 GMT. 
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Table X. Atmospheric Parameters around White Sands. 


Atmospheric Parameter 

Value 

Units 

Ambient Air Temperature 

35 

°G 

Surface Air Temperature 

75 

°G 

Ambient Relative Humidity 

.33 

fraction 

Surface Relative Humidity 

1.0 

fraction 

Ground Type 

sand 

~ 

Wind Speed 

3.1 

m/s 

Wind Direction 

27° 

NNE 

Surface Pressure 

1003 

hPa 

Visibility 

16 

km 

Lower Level Glouds 

2000 

m 

Lower Level Gloud Type 

Scat 

-- 

Mid Level Glouds 

none 


Upper Level Glouds 

3700 

m 

Upper Level Gloud Type 

Scat 

- 


C. SELECTION OF TARGET LOCATIONS 

In each case studied, the target location entered into TAWS corresponded to 
a coordinates of a nearby region. These coordinates are summarized in Table XII. 


D. PROCEDURES 

In this study, the MRTD was calculated for chosen atmospheric conditions 
using both FLIR92 and NVThermIP. 

1. Parameters Used in Calculations 

The parameters used for the FLIR92 and NVThermIP calculations are listed 
in tables in Appendix D. Whenever possible, the same inputs were used in FLIR92 
and NVThermIP. 

2. FLIR92 Calculations 

The FLIR92 calculations were two-fold. First, a MatLab version of the cur¬ 
rent operational FLIR92, supplied by Dr. Andy Goroch, was used. Then the MatLab 
version of FLIR92 was modihed to include atmospheric turbulence. In order to deter- 


60 




Table XL Atmospheric Parameters around Point Conception. 


Atmospheric Parameter 

Value 

Units 

Ambient Air Temperature 

25 

°C 

Surface Air Temperature 

31 

°C 

Ambient Relative Humidity 

.55 

fraction 

Surface Relative Humidity 

1.0 

fraction 

Ground Type 

water 

- 

Wind Speed 

3.7 

m/s 

Wind Direction 

290° 

WNW 

Surface Pressure 

1016 

hPa 

Visibility 

16 

km 

Lower Level Clouds 

none 


Mid Level Clouds 

none 

- 

Upper Level Clouds 

none 

- 


Table XII. Target Locations. 

The target coordinates are entered into TAWS. 


Location 

Coordinates 

N White Sands Range 

33°47T0”N;106°27’47”W 

Bnoy 46063, Pt Conception 

34°16’21”N;120°4F55”W 


mine the Monin-Obukhov length and T*, g*, and n*, certain atmospheric conditions 
were entered into the MatLab FLIR92 program, as ontlined in Table XIII below. 
Most actnal valnes depended on the atmospheric conditions and are listed later for 
the specihc cases. 


Table XIII. Atmospheric Condition Variables for Input into Modified FLIR92. 


Variable 

Units 

Air Temp 

°C 

Sfc Temp 

°c 

Rel Humidity Air 

fraction 

Rel Humidity Sfc 

fraction 

Wind Speed 

m/s 

Sfc Pressnre 

hPa 

Thermal Sfc Roughness 

0.2 m 

Momentum Sfc Roughness 

0.2 m 
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3. NVThermlP Calculations 

No modifications to NVThermlP were made, but the output was modihed as 
discussed in Section I.l.C. 

4. TAWS Calculations 

The resultant MRTD curves were entered into TAWS version 5.2 and the 
maximum range detection results from TAWS were compared. The target used in 
TAWS was the 155 mm M795, a projectile in a black plastic bag. The sensor used in 
TAWS was the generic 1000 IR sensor, modihed with the MRTD outputs from either 
FLIR92 or NVThermlP. Weather inputs are given in Tables X and XL The test cases 
were run at a height of 100 m. 

For the MRTD input, a table of spatial frequency and corresponding MRTD 
result was inserted into TAWS. The FLIR92 output had a larger frequency interval 
than NVThermlP. To account for this, the NVThermlP MRTD results were inter¬ 
polated onto the FLIR92 results. Thus, the same number of frequency-MRTD pairs 
were entered for both the NVThermlP and FLIR92 cases. 
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V. 


RESULTS 


A. ATMOSPHERIC TURBULENCE IN FLIR92 

The modified FLIR92 including an MTF for atmospheric turbulence with a 
varying dependence on was compared to the original FLIR92. The inputs are 
given in Appendix D with the atmospheric inputs for the modified FLIR92 given in 
Section IV.B for White Sands, NM. The resultant MRTDs and system MTFs were 
compared in the horizontal and vertical directions, as shown in Tables XIV and XV. 
M indicates the modified FLIR92 with atmospheric turbulence, O the original, fs is 
in mrad“^. 

Table XIV. MRTD for Modified, Original FLIR92 in the horizontal direction. 
The M indicates the modified FLIR92 while the O indicates the original FLIR92. fg 
is in mrad“^ and the MRTD is in K. 


fs 

MRTDm 

MTFgys^ 

MRTDo 

MTFgys^ 

0.9935 

0.0013 

0.8214 

0.0013 

0.8214 

2.9804 

0.0054 

0.4874 

0.0054 

0.4875 

4.9674 

0.0155 

0.2429 

0.0155 

0.2429 

6.9543 

0.0453 

0.1008 

0.0453 

0.1008 

8.9413 

0.1511 

0.0343 

0.151 

0.0343 

10.9283 

0.6072 

0.0094 

0.6069 

0.0094 

12.9152 

3.13 

0.0019 

3.1284 

0.002 

14.9022 

23.2345 

0.0003 

23.2193 

0.0003 


Figure 19 illustrates how the difference between the modified and original 
FLIR92 is only apparent at higher frequencies. Even so, the scale indicates how small 
is the impact of including atmospheric turbulence in the MRTD calculation. 

There is a difference in resultant MRTD when including the atmospheric tur¬ 
bulence MTF, as seen at the higher frequencies, so starting at 10.9283 cycles/mrad 
in the horizontal case. The difference is minimal, though, less than 0.2% of the hor¬ 
izontal MRTD and vertical MRTD. The difference is too small to impact the range 
calculation in TAWS since the input for MRTD in TAWS is only of accuracy 10“^ K. 
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Table XV. MRTD for Modified, Original FLIR92 in the vertical direction. 
The M indicates the modihed FLIR92 while the O indicates the original FLIR92. 
is in mrad“^ and the MRTD is in K. 


fs 

MRTD]\^ 


MRTDo 

MTF^yso 

0.9935 

0.0013 

0.8254 

0.0013 

0.8254 

2.9804 

0.0053 

0.5094 

0.0053 

0.5095 

4.9674 

0.014 

0.2751 

0.014 

0.2751 

6.9543 

0.0365 

0.1295 

0.0364 

0.1295 

8.9413 

0.1023 

0.0528 

0.1023 

0.0528 

10.9283 

0.3224 

0.0185 

0.3222 

0.0185 

12.9152 

1.1786 

0.0054 

1.178 

0.0054 

14.9022 

5.2264 

0.0013 

5.223 

0.0013 

16.8891 

30.9544 

0.0002 

30.9296 

0.0002 


Again, Figure 20 indicates the differences at higher frequencies and the small scale of 
the differences. 


Horiz: Modified FLiR92 - Originai FLIR92 MRTD 



Figure 19. Horiz Difference Between MRTD Curves for Modihed, Original FLIR92. 
Differences are not observed until higher frequencies. 
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B. COMPARING FLIR92 AND NVTHERMIP 
1. Comparing Resultant MRTDs 

The modified FLIR92 MRTD, including an MTF for atmospheric turbulence 
with a varying dependence on was compared to NVThermIP. As explained in 
Section I.l.C, the CTF results from NVThermIP were used to find a comparable 
MRTD. From both FLIR92 and NVThermIP, the MRTD results were split into hori¬ 
zontal and vertical MRTDs. Comparisons of the horizontal MRTD show fairly similar 
results, as seen in Figure 21. The downturn on the FLIR92 MRTD curve is due to 
insufficient eye filtering. NVThermIP attempts to compensate for this weakness, and 
the difference is apparent in the upturning curve at low frequencies. Overall, the 
NVThermIP MRTD is greater because NVThermIP includes an improved and less 
pessimistic calculation to incorporate the human visual response to the display (see 
Section B.l.n). 

Comparing the vertical MRTD, as in Figure 22, shows that once again at low 


Vert: Modified FLiR92 - Originai FLiR92 MRTD 



Figure 20. Vert Difference Between MRTD Curves for Modified, Original FLIR92. 
Differences are not observed until higher frequencies. 
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spatial frequencies, FLIR92 turns down while NVThermlP turns up. At high spatial 
frequencies, though, the two curves appear to agree closely in terms of slope and 
value. Because this study focuses on differences between the two MRTD models and 
the horizontal MRTD exhibits a greater difference, the horizontal MRTD is used as 
the input into TAWS. 

2. Comparing Detection Range 

The MRTDs were calculated with FLIR92 and NVThermlP for both White 
Sands and Point Conception. 

The resultant maximum range detection from TAWS at 35%, 65%, and 95% 
probabilities, are compared for both the White Sands site and the Point Conception 
site. NVThermlP uses a constant atmospheric structure parameter. was set 



Figure 21. Horizontal MRTD Curves for FLIR92, NVThermlP. 

The horizontal MRTD curves have different behaviors at the low frequencies and 

high frequencies. 
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to m“i for both cases. Other inputs to TAWS are listed in Section IV.B. 

The White Sands case shows NVThernilP maximum detection range going from 0% 
to 1.5% greater than the FLIR92 maximum detection range, for all probabilities. 
In Figure 23, the maximum detection ranges at all probabilities for MRTD from 
NVThermlP (top) and FLIR92 (bottom) are shown. 

Figure 24 better illustrates the differences between the maximum detection 
ranges calculated using the MRTDs from NVThermlP and FLIR92. The largest 
percent difference is at 35% probability of detection. 

The Point Conception case shows NVThermlP maximum detection range go¬ 
ing from 0% to 2% greater than the FLIR92 maximum detection range, for all prob¬ 
abilities. Figure 25 illustrates the maximum detection ranges for 35%, 65%, and 95% 
probabilities of detection based on the MRTDs from both NVThermlP (top) and 
FLIR92 (bottom). 



Figure 22. Vertical MRTD Curves for FLIR92, NVThermlP. 
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To better examine how the maximum detection ranges based on MRTDs from 
FLIR92 and NVThermlP compare, Figure 26 shows FLIR29 and NVThermlP max¬ 
imum detection range results for each probability. The largest percent difference is 
again at 35% probability of detection. 

3. Comparing Marine/Desert Conditions 

After comparing maximum detection ranges based on MRTD calculations from 
FLIR92 and NVThermlP, differences between maximum detection ranges at White 
Sands and Point Conception were compared. As shown in Figure 27, at 35% prob¬ 
ability of detection, there is a greater difference between maximum detection range 
from the FLIR92 MRTD and the NVThermlP MRTD for Pt Conception. At 95% 
probability of detection, as shown in Figure 27, there is much less difference between 
maximum detection ranges using FLIR92 and NVThermlP MRTDs, but there is still 
an observeable difference. 

The difference between NVThermlP and FLIR92 at each location is better 
seen in the bar plots in Figures 28 and 29 showing the viewing angle versus the 
maximum detection range at 35% and 95% probabilities of detection respectively. 
The 35% case most clearly demonstrates that at White Sands, over the hot desert, 
there is the greatest difference in the maximum detection range between FLIR92 
and NVThermlP, with NVThermlP having a greater range. At 95% probability of 
detection, the difference between the desert and the marine locations is much smaller. 




Figure 23. Comparing Maximum Detection Range at White Sands using FLIR92, 
NVThermIP. 

Polar Plot of viewing angle (deg) versus maximum detection range (km) for 35%, 

65%, and 95% probabilities of detection, using the MRTD from NVThermIP (top) 

and FLIR92 (bottom). 
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Max Det Rng 35%, 65%, 95% Prob White Sands 


WS Comparing Max Detection 35% Prob 


Comparing Max Detection 65% Prob 
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Comparing Max Detection 95% Prob 
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Figure 24. Maximum Detection Range for White Sands at Different Probabilities. 
Polar Plot of viewing angle (deg) and maximum detection range (km) for both 


FLIR92 and NVThermlP. 
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Figure 25. Comparing Maximum Detection Range at Pt Conception using FLIR92, 
NVThermIP. 

Polar Plot showing view angle (deg) versus the maximum detection range (km) for 
35%, 65%, and 95% probabilities of detection using MRTD from NVThermIP (top) 

and FLIR92 (bottom). 
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Max Det Rng 35%, 65%, 95% Prob, Pt Conception 


PC Comparing Max Detection 35% Prob 

180 c 



Comparing Max Detection 65% Prob 
180 , 



Comparing Max Detection 95% Prob 
180 4 




Figure 26. Maximum Detection Range for Pt Conception at Different Probabilities. 
Polar Plot showing view angle (deg) versus the maximum detection range (km) for 


both FLIR92 and NVThermIP. 
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Figure 27. Comparing Maximum Detection Ranges at White Sands and Pt Concep¬ 
tion at 35% Prob and 95% Prob. 

Polar Plot showing view angle (deg) versus the maximum detection range (km) for 
both FLIR92 and NVThermlP at both White Sands (WS) and Pt Conception (PC) 
at 35% and 95% probabilities of detection. 


73 



























5 


White Sands: NVThermlP and FLIR92 Max Det Rng 35% 



Pt Conception: NVThermlP and FLIR92 Max Det Rng 35% 



Figure 28. Comparing Maximum Detection Ranges at White Sands and Pt Concep¬ 
tion at 35% Prob. 

Bar plot showing view angle (deg) versus the maximum detection range (km) for 
both FLIR92 and NVThermlP at White Sands and Pt Conception at 35% 

probability of detection. 
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White Sands: NVThermlP and FLIR92 Max Det Rng 95% 



Viewing Angle (deg) 


Pt Conception: NVThermlP and FLIR92 Max Det Rng 95% 



Viewing Angle (deg) 


Figure 29. Comparing Maximum Detection Ranges at White Sands and Point Con¬ 
ception at 95% Prob. 

Bar plot showing view angle (deg) versus the maximum detection range (km) for 
both FLIR92 and NVThermlP at White Sands and Pt Conception at 95% 

probability of detection. 
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VI. CONCLUSIONS AND 
RECOMMENDATIONS 

This study considered two cases to test possible improvements to the MRTD 
input into TAWS. First, FLIR92 was modified to include atmospheric turbulence and 
the resultant MRTDs from the original and modified FLIR92 were compared. Initial 
results indicate that while there is a difference in including atmospheric turbulence 
in the MRTD calculation, the difference in most cases is too small to be significant. 
Second, FLIR92 and NVThermlP MRTD results were entered into TAWS and the 
resultant maximum detection ranges were compared. Again, the differences in the 
representative test cases were minimal. In comparing the FLIR92 and NVThermlP 
results in TAWS, the percentage improvement in using NVThermlP was slightly 
greater for the marine test case than for the desert test case. 

A. CONCLUSIONS FLIR92 WITH/WITHOUT ATMO¬ 
SPHERIC TURBULENCE 

The results using a test case as described in Chapter IV indicates very little 
difference in the MRTD when an atmospheric turbulence MTF is included in the total 
system MTF. The atmospheric turbulence MTF is very close to unity, on the order 
of 1 — 10“^, so the total system MTF is simply being multiplied by a factor of unity. 
In the Goodman long exposure MTF given in Equation A.103, the general form is an 
exponential that depends on fg, A, and C^{z) where is the height in the atmosphere. 
Three possible cases exist where the atmospheric MTF may become significant. First, 
if the spatial frequency is very large, corresponding to imaging very small targets, the 
MTF may be significant. Second, if the diffraction wavelength is very small, well 
below the bounds of an IR study, the MTF may be significant. Finally, if is large, 
at a minimum on the order of 10“^'^ m“i, then the atmospheric MTF may become 
more important. 
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Near the surface where the bulk model is applicable, depends on the am¬ 
bient air and surface temperatures, wind speed, surface pressure, relative humidity, 

4 

thermal surface roughness and momentum surface roughness as well as z~3. Figure 
30 shows that for the White Sands case near the surface, is slightly greater than 
IQ-io indicating that only very close to the surface may the atmospheric MTF 
be important. Too close to the surface, though, at less than 0.2 m, the bulk model 
for including approximations for T*, u*, and g* may no longer be valid due to the 
extreme heat fluxes and bumps along the surface itself that are likely greater than 
the sensing height and will affect the scaling parameters. 


Bulk Model Near Sfc T , = 35°C, T , = 75°C 
amb sfc 



Figure 30. Near Surface Behavior of the Bulk Model. 
The height above the surface shown here is 0 to 1 pm. 


B. CONCLUSIONS FLIR92 VERSUS NVTHERMIP 

Comparing the maximum detection range using the horizontal MRTD from 
FLIR92 and NVThermlP, there is not a signihcant difference between the results 
in the marine environment and the results in a desert environment. In both cases. 
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the MRTD from NVThermlP yielded a greater maximum detection range by up to 
2%, but the results did not even consistently improve by 2%. Figures 21 and 22 
show substantial difference between FLIR92 and NVThermlP MRTD results, but 
apparently these differences do not signihcantly impact TAWS calculations. 

C. RECOMMENDATIONS 

Since neither the modification of FLIR92 nor the use of NVThermlP instead 
of FLIR92 resulted in signihcant differences in TAWS, it seems best at the moment 
to continue using FLIR92 to calculate the MRTD for TAWS. However, should TAWS 
increase its sensitivity to the MRTD, the differences described in this study may 
become signihcant. 

Including an atmospheric turbulence MTF in FLIR92 simply involves multi¬ 
plying the total system MTF by the atmospheric MTF. Such a simple adjustment 
may prove useful for cases where is large. The same adjustment could be made 
in NVThermlP, where the current atmospheric MTF assumes a constant C"^. In NV¬ 
ThermlP, the maximum is m“i but using such a large up through the 

atmosphere is highly unrealistic. 

If the use of FLIR92 for MRTD calculations for TAWS continues, other mod- 
ihcations will need to be made. FLIR92 assumes a CRT display, but in many ap¬ 
plications now the images are displayed on a hat panel. It is possible to enter a 
user-dehned display MTF describing the hat panel, but updating the program would 
be a better solution. 
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APPENDIX A. GOODMAN DERIVATION FOR 

LONG EXPOSURE 

The following derivation is primarily summarized from Goodman (1985), with 
substantiations from other sources. 

A generally accepted form of the refractive index of air at optical frequencies 
is given by Goodman as 

p 

n = l + 77.6(1 + 7.52 x lO’^A-^)- x 10"® (A.67) 

where A is the wavelength of light in /rm, P is the atmospheric pressure in mb, and 
T is the atmospheric temperature in K. To simplify, consider that pressure variations 
are small compared to the temperature fluctuations at optical frequencies and so may 
be neglected. Then, choosing a representative A = 0.5/im, it is apparent that for an 
incremental change dT of temperature, 

79P 

dn = —— X 10~^dT. (A.68) 

Inhomogeneities in the refractive index, otherwise described as turbulent ed¬ 
dies, can be considered individual packets of air with individual characteristic refrac¬ 
tive indices. Then, the power spectral density of homogeneous turbulence, <hn(^), 
where k, is the vector wavenumber k, = {kx, i^z), yields a measure of the relative 
number of eddies of size Lx = — and Lv = — and Lz = —■ For isotropic turbu- 
lence, $„(«;) is only a function of the non-vector wavenumber k that corresponds to 
an eddy of size L where L = —. 

In the classic theory of turbulence outlined by Kolmogorov (1961) and ex¬ 
panded by Goodman, the power spectral density consistss of three different regions 
depending on the size of k. Kolmogorov hypothesizes that the kinetic energy associ¬ 
ated with the larger eddies will be redistributed to successively smaller and smaller 
eddies without loss of energy. For the smallest k, corresponding to large scales since 
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K, = 1/A, turbulence depends on large-scale topographic and meteorological condi¬ 
tions. This makes it unlikely that the turbulence would be isotropic or homogeneous, 
so no mathematical form of turbulence is predicted in this region. Once k grows 
beyond some critical wave number kq, *h„(/T;) can be determined by physical laws 
governing how the large turbulent eddies break up into smaller ones. In this re¬ 
gion beyond Kq, the so-called inertial subrange of the spectrum, the power density 
spectrum may be written as 

(A.69) 

where again C\ is the atmospheric structure constant. Goodman further dehnes 
C\ as the structure constant of the refractive index fluctuations and a measure of 
the strength of the fluctuations. In the third region, beyond another critical wave 
number Km-, turbulent eddies are very small scale and viscous forces may cause them 
to dissipate their energy. This inner scale of the turbulence results in rapid decay of 
where n > nm- The decay can be modeled as 


^n{K) = 3 exp (A-70) 

Tatarski (1961) points out that this model is only valid when for the inner scale of 
the turbulence, /q = ^, is chosen equal to and k> Kq. The above expression 
approaches inhnity at the origin, so Goodman uses the von Karman spectrum instead. 
This form also accounts for the hnite amount of air in the atmosphere and does not 
become arbitrarily large as k goes to 0. Approximated, the von Karman spectrum 
yields a power spectrum density of 



30(^2 -|- Kq 


(A.71) 


and Goodman emphasizes that this mathematical form is only an artihcial means of 
avoiding the poles. The von Karman spectrum depends on the critical wave nnmbers 
in both the medium and small size turbulent eddies. Figure 31 illustrates the form 
of the von Karman spectrum compared to the Kolmogorov spectrum. 
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Kolmogorov vs von Karman Spectrum 



Figure 31. The Kolmogorov and von Karman Spectra. 

The Kolmogorov spectrum approaches inhnity as k goes to zero, but the von 
Karman spectrum avoids the division by zero and has a hnite value at k = 0. 


The structure function of the refractive index fluctuations in the case of at¬ 
mospheric turbulence in imaging systems is dehned as 

Dn{fl, r^) =< [ni(rl) - ni(r^)]^ > (A.72) 

where ni is the random refractive index and r = rl — = (Ax, Ay, Az). Assuming 

homogeneous turbulence, 

Ii„(f) = 2[r„(0)-r„(f)] (A.73) 

where the autocorrelation function T^ may be expressed as 

OO 

rn(r) = JJJ exp {-jK ■ f)d^K (A.74) 

— OO 

and it follows that 

OO 

r„(0)= JJJ (A.75) 
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Combining, then, 


DJf) = 2 


[1 — exp {—jn ■ r)] 


(A.76) 


Assuming a real solution so the exponential simplihes to cos{K-r), and when r„(r') is 
spherically symmetric and ni is statistically isotropic, the Fourier transforms of the 
power density function and the autocorrelation function may be expressed as single 
integrals. In this case, for k = ^+ and r = + {5yY + {Sz^, 

A °° 

Att 


47r f 

Trifh) = — / sin(fi;r)(i/T;. 

r J 
0 

Then the structure function simplihes: 

OO 

Dn{r) = 2r„(0) — Stt [ ^nM —sin(Kr)dK. 

J r 

0 

To determining r„( 0 ), the limit as r —0 is taken. 


, ^ / N . . /■ ^ / X oSmiKr) , 

hmr„(r) = lirndvr / - dn. 


(A. 77) 


(A.78) 


r —>0 


r^O 


nr 


(A.79) 


Then, since the integral does not depend on r, the limit may be moved inside the 
integral to hnd 


limTnfr) =471 <F„ ( k) lim 


sm /tr 


1—>0 


r^o nr 


dn. 


Taking the limit, as r approaches 0 is found to be 

OO 

limT^fr) =471 
r^O J 

0 

which yields a Dn{r) of 


(A.80) 


(A.81) 


Dn{r) = Stt y $„(«;) 
0 


1 - 


sm Kr 


nr 


dn. 


(A.82) 


Assuming that the power density has the form of Equation A.69, the structure func¬ 
tion may be simplihed to 




1 - 


sm Kr 


nr 


dn. 


(A.83) 
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Equation A.83 is integrated using the integral identity 


X 


sin tax j 
ax 


dx = 


T(z/) sin( 


TTU ^ 
2 ^ 


^1+u 


where —3 < v < —1. It is valid to use Equation A.84 because in this case v 


D (r) = — 

nl ) gu 


5 


r ( ) sin ( j 


SttX 




D„{r) = lx 


(A.84) 


= -!> so 

(A.85a) 

(A.85b) 


where the coefficient of unity in Equation A.85b is due to the chosen dehnition of 

Goodman’s derivation continues as follows by determining the optical transfer 
function (OTF). The derivation will be summarized here, but the specihc details are 
available in Goodman’s text (1985). The MTF for atmospheric turbulence is simply 
the modulus of the OTF, and can easily be calculated. In this case, the OTF is a 
real value, so the modulus of the OTF is the same as the OTF, so the MTF equals 
the OTF. This derivation focuses on long-exposure systems where the exposure time 
is much greater than one one-hundredth of a second. The short-exposure case can 
also be signihcant in applications but is neglected in this study. Also, the long- 
exposure OTF is assumed to be identical to the ensemble-averaged OTF. Using this 
assumption, Goodman derives a total average OTF of the form 


OTF{p) = OTFo{p)OTFl{v) 


(A.86) 


where v = y and vu and z/y are spatial frequency variables. To simplify the 

derivation, a circularly symmetric OTF is assumed, and the long-exposure OTF of 
the atmosphere, OTF^, has the form. 


OTFl = exp 




(A.87) 


where again D{\fu) is the wave structure function, A = | where c is the speed of 
light, and 9 is the center frequency of the waveband under consideration. 
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In the near-field calculation of the wave structure function, a detailed expres¬ 
sion for the wave structure function is necessary to specifically define the form of the 
long-exposure atmospheric OTF. Three assumptions are made. 

First, the target is assumed to be very far away from the detector lens. In 
this way, the angular extent of the target is small enough that over the long time 
average, the entire target is affected uniformly by the atmosphere. This assumption is 
considered the ’isoplanatic’ assumption, in other words assuming that the point spread 
function is uniform over the entire FOV. In the long-time-average, this assumption is 
probably a close approximation to truth and thus is not very restrictive. 

Second, the atmospheric turbulence is assumed to be homogeneous and isotropic 
over a finite distance in front of the detector imaging lens. This assumption is in¬ 
valid when considering vertical viewing, which is a problem for situations such as 
an aircraft imaging a ground-based target. This assumption will be removed later 
though, as varies over the pathlength. 

Third, the imaging system is assumed to be so far within the near field of 
the relevant turbulent eddies that each ray incident on the inhomogeneous medium 
may be approximated to merely be delayed by that medium with no bending of the 
rays. In other words, the turbulence is considered weak enough that no significant 

p 

amplitude scintillation effects occur. This assumption is only valid when z ^ 
Clearly, in actual practice this assumption is also unrealistic, but will be removed 
later. 

Without going into the details of the derivation, these assumptions are used 
to find that the wave structure function equals the phase structure function and 
the atmospheric OTF may be equivalently expressed in Equation A.87 and below in 
Equation A.88, 

OTFl = exp [-1 d5(A/z/)] (A.88) 
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where of course Ds{\fi^) is the phase structure function. Using the covariance func¬ 
tions Cn, not to be confused with the atmospheric structure constant C^, 

z oo 

Ds{r) = ¥ J dz' J 2Cco {z' - z") - 2Cn\l{z' - z"f + dz" (A.89) 

0 0 

where again is the distance from the target to the imaging lens. The difference of the 
covariance functions may further be expressed as a difference of structure functions, 
Dn, 

Z Z _ 

Ds{r) = F J dz' J Dn\/ {z' - z'f + {z' - z") dz". (A.90) 

0 0 

Using the simplihcation that the integrand is an even function of z' — z", the double 
integral may be reduced to a single integral such that with = z' — z", 

Z 

Ds{r) = 2F J {z — Az) Dn\J (Az)^ -|- (Az) — d{Az). (A.91) 

0 

Next, the specihc form of the structure function that was found above. Equa¬ 
tion A.85b, is substituted into Equation A.91 to yield an expression for the phase 
structure function, 

Z 

Ds{r)=2FCl j{z-Az)[{[Azf+ -[Az^]d{Az). (A.92) 

0 

Because Equation A.85b is only valid for Iq < r < Lq, Equation A.92 is only valid for 
A^ < Lq. This makes it appear that Equation A.92 is only valid when the pathlength 
is less than the outer scale Lq, but in the typical cases where the path separation 
r is less than Lq, the integrand vanishes for large Az and it does not matter what 
was the exact form of the structure function for r > Lq. For Az > Lq and r Lq, 
the integrand approaches zero. Thus, Ds{r) can be used for all pathlengths without 
violating any already stated assumptions. 

Since Az is much less than z. Equation A.92 may be rewritten, neglecting the 
term of Az multiplied by the rest of the integrand, 

OO 

Ds{r) ^ 2FcIz j [{[Az]^ + FF " [Az]i]d(Az) (A.93) 

0 
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where the upper integration bound has been extended to cx) since from to cx) barely 
contributes to the integrand. This integral may be evaluated numerically using a pro¬ 
gram such as Mathematica. First, a change of variables makes the integral expression 
clearer: Az = ru and dAz = rdu, 

OO 

Ds{r) = Ds{r) = PC^z J {P + Ip — du. (A.94) 

— OO 

Mathematica gives the result that 

Ds{r) = 2.9U3PClzr^. (A.95) 


Substituting Equation A.95 into the expression for the long-exposure atmospheric 
OTF, Equation A.88, the hnal expression for the long exposure atmospheric OTF is, 

OTFi(z/)=exp —^x2.9lPC'^z(\fh'"j^ . (A.96) 

To obtain a form of OTFl that is independent of the optical system, the OTF 
is expressed in terms of the frequency fg measured in cycles/rad rather than u in 
cycles/m. In that case, fs = fv and 

OTFPfg) = exp [-i X 2.9lPClz (A/,)"] . (A.97) 

To further simplify, k = ~, so Equation A.97 becomes 

OTFPfg)=exp -5.82 ^ 

As 

The long exposure OTF is identical to the long exposure MTF since the MTF is 
merely the modulus of the OTF and in this case the modulus is equal to the OTF 
itself. Also, though the above result is limited to the three assumptions stated above, 
it may be generalized. 

First, the assumption that refractive index perturbations only delayed rays 
rather than bending or diffracting them may be eliminated without any change to 
the long exposure OTF theory. Goodman proceeds through an analysis considering 
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(A.98) 



the case of refractive index perturbations occurring within a finite region along the 
propagation path and then integrates along the entire propagation path. This gener¬ 
alization leads to exactly the same results as for the previous analysis, indicating that 
it is possible to consider cases where the rays diffract and bend as they pass through 
varying refractive indices. 

Second, since variations in occur frequently in practice, the assumption that 
the turbulence is uniform is removed. The analysis in Goodman examines smooth 
and slow variations to the atmospheric structure constant over the path length in 
both the vertical and horizontal directions. Goodman points out that his analysis is 
an approximation, but has been verihed with more rigorous methods. Since Good¬ 
man assumes a slow and smooth change in G^, his method of determining the OTF 
considers a ” quasihomogeneous” model for the structure function of the refractive 
index fluctuations, 

G'n(Gi,r^) = Cl \n - G 2 I 3 (A.99) 

where /q < |G ~ gI < -^o and G and are vectors indicating the location of the 
changes in G^. Glearly, in this assumed structure function of the refractive index 
fluctuations, signihcant changes in Cl only occur over distances greater than or equal 
to Lq. Assuming that the atmosphere can be divided into stacking layers of thickness 
Az along the propagation path such that the fluctuations in log amplitude and phase 
in the different layers are uncorrelated, the wave structure function after passing 
through N layers may be expressed as, 

N 

D{r) = Y.D,{r) (A.IOO) 

i=l 

where Di is associated with the Ah layer. Taking Zi to represent the 2 : coordinate in 
the middle of the Ah layer and using Equation A.97 for each individual layer wave 
structure function, the total wave structure function may be expressed as 

N 

D{r) = 2.9lP^G2(;2i)Azri. 

i=l 
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Generalizing because variations of are assumed to be slow compared to the length 
A;^, replace the finite sum with the integral along the propagation path, 

Z 

D{r) = 2.9lPrt J Cl{z')dz' (A.102) 

0 

where z' is the variable along the propagation path from 0 to 2 :. Using Equation A.102 
in the expression for the long exposure OTF, Equation A.88, the long-exposure OTF 
with a varying is 

Jcl{z')dz' 

OTFz(f,)=exp -5.82irV/S-,- . (A.103) 

As 

See Section III.A for details on the approximations for how G^ varies in the 
vertical. 
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APPENDIX B. EXACT FORM EYE/BRAIN 
SPATIAL INTEGRATION 


The complete expressions for the eye/brain system spatial integrations are 
given below, assuming the eye/brain spatially integrates over all frequencies within 
the yet undegraded image of the target. The hrst subscript in these equations is the 
direction of spatial integration. The second is the orientation of the MRTD target. 

In scanning systems, for horizontal MRTD target patterns, the system spatial 
integrals are given in the horizontal and vertical spatial integration directions. In the 
horizontal. 


EkAf.) 


JP Q 



2 /. 


(B.l) 


and in the vertical. 



— OO 



2fs 


(B.2) 


where and are the horizontal and vertical system noise hlter MTF 

components as given in Equation 2.20, A/p is dehned by Equation 2.24, is the 
samples per detector vertical IFOV and is the detector vertical IFOV in mrad. 

For vertical MRTD target patterns in scanning systems, the system spatial 
integrals are given in the horizontal and vertical spatial integration directions. In the 
horizontal. 


EhMs) 


Vs 

^fp 


OO 

1 

B' 

_ 1 

Ttzi/ 

J 

0 

Vs 


dv 


(B. 3 ) 


and in the vertical. 



— OO 



Vs 


(B.4) 
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For horizontal MRTD target patterns in staring systems, the system spatial 
integrals are given in the horizontal and vertical spatial integration directions. In the 
horizontal, 


EkU.) 



(B.5) 


and in the vertical. 



— OO 



2/« 


(B.6) 


where 6h is the detector horizontal IFOV in mrad and Sh is the samples per detector 
horizontal IFOV. 

For vertical MRTD target patterns in staring systems, the system spatial in¬ 
tegrals are given in the horizontal and vertical spatial integration directions. In the 
horizontal. 


EkU.) 




2 /. 


(B.7) 


and in the vertical. 



— OO 



2fs 


(B.8) 
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APPENDIX C. TEMPERATURE SCALING OF 

O'TVH 

In order to scale (Jtvh, the model scales the detector spectral detectivity from 
300 K to Ti where Ti is the new temperature. This scaling can only occur if the 
detector noise is limited by the incident background photon flux, in other words, the 
background limited in performance (BLIP) condition. The BLIP conditions happens 
when the background photon flux is much greater than the signal flux, causing photon 
noise to be the dominant noise factor (Driggers 1999). The scaling of the detector 
spectral detectivity is given in Equation C.l, 

OO 

/QA(A,300)dA 

D*{X, Ti) = D*{X, 300) - (C.l) 

. /QA(A,Ti)dA 
^ 0 

where Qv{X,T) the spectral photon flux density where T = 300,Ti. Q^{X,T) is 
defined in photons/s/cm^. 

For BLIP detectors, using this scaled D*{X,Ti), and for non-BLIP detectors 
using the unsealed D*, the appropriately temperature-scaled atv is given as 

/D-(A,300)Jg5(A)dA 

{o'tvh)Ti = {O'tvh)300^ -• (C-2) 

fD*{X,T,)-^{X)dX 


93 



THIS PAGE INTENTIONALLY LEFT BLANK 


94 



APPENDIX D. TABLES OF INPUTS TO 
FLIR92, NVTHERMIP 


The salient inputs to NVThermlP are listed in Tables XVI, XVII, and XVIII. 
Inputs to NVThermlP are made in the form of a text hie, so in the text hie each entry 
is of the form ‘variable,value’ as in ‘ SensorName, Sensor Test 1’ if SensorName is the 
variable and SensorTestl is the value, or the name of the sensor being tested. Inputs 
to FLIR92 are listed in Tables XIX and XX. These inputs are made in a MatLab hie. 
The values given are generic and correspond to no special case. The H stands for the 
horizontal direction and the V for the vertical direction. Atmospheric data changed 
for each case and is listed in a separate Tables X and XI given in Chapter IV. Note 
that calculations were carried out for a staring sensor. 


Table XVI. Input Parameters NVThermlP Part I. 


Variable 

Value 

Units 

Comments 

SensorName 

SensorTestl 

- 

Any name okay 

TypeOfImager 

0 

- 

0 for Staring 

SingleFrame 

1 

- 

- 

Cut OnWavelength 

3.3 

fj,m 

- 

Cut OffWavelength 

5 

fj,m 

- 

Magnihcationin 

0 

- 

0, Calculated in NVThermlP 

HorizontalFOV 

0.48 

deg 

Same as FOV calc, in FLIR92 

VerticalFOV 

0.48 

deg 

Same as FOV calc, in FLIR92 

FrameRate 

30 

Hz 


Verticalinterlace 

2 

- 

- 

Electronicinterlace 

1 

- 


DiffWavelength 

4.15 

/im 

~ 

AveOpticTrans 

1 

frac 

- 

FocalLengthIn 

45.7 

cm 

- 

FNumberIn 

0 

- 

0, Calculated in NVThermlP 

Apert ureDiameterIn 

15.2 

cm 

- 

OpticsBlur 

0.01 

mrad 

- 

OpticsBlurUnits 

0 

- 

0 for object space 

OpticsBlurType 

0 

-- 

0 for RMS 

VibBlurX 

0.02 

mrad 
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Table XVII. Input 


Variable 

Value 

VibBlurY 

0.02 

VibBlurType 

0 

DetHorrDim 

23 

DetVertDim 

23 

PeakDstarIn 

lE+12 

IntegrationTimeIn 

7000 

NumTDI 

1 

ScanEfficiency 

1 

NumHorrDet 

512 

NumVertDet 

512 

SamplesPerHIFOV 

1.4 

NoiseFactorX 

0 

NoiseFactorY 

0 

Sigmavh 

0 

Sigmav 

0 

Sigmah 

0 

FPN 

0 

PtSi 

1 

EmissionCoef 

0 

BarrierHeight 

0 

Uncooled 

1 

Dither 

0 

LowPassSdB Cutoff 

22000 

LowPassF ilter Order 

1 

NoiseGain 

0 

Frameintegration 

1 

Interph 

0 

Interpv 

0 

InterpTypeh 

0 

InterpTypev 

0 

EZoom 

0 

EZoomType 

0 

Boosth 

0 

Boostv 

0 

DispType 

0 

EoMUX 

1 

HorLEDSize 

0 

Vert LED Size 

0 


Parameters NVThermlP Part II. 


Units Comments 


mrad 



0, Calculated 


0, Calculated 


0, Calculated 


0, Calculated 


0, Calculated 


0 , none 


1 , no 


0, Calculated 


0, Calculated 


1 , no 


0 , no 


in NVThermlP 


in NVThermlP 


in NVThermlP 


in NVThermlP 


in NVThermlP 


in NVThermlP 


in NVThermlP 


Number hlter poles 



0, no H interopolation 


0, no V interpolation 


0, no H interpolation 


0, no V interpolation 


0, no EZoom 


0, no EZoom 


no H boost 


no V boost 


0, CRT 


1 , no 
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Table XVIII. Input Parameters NVThermlP Part III. 


Variable 

Value 

CRTType 

0 

LEDHeight 

0 

LEDWidth 

0 

DispSpotHeight 

0.02 

DispSpotWidth 

0.02 

AveDispLum 

10 

MinDispLum 

0 

SceneContTemp 

10 

DispHeight 

15.24 

DispViewDist 

70 

NumEyes 

2 

At mTr ansTyp e 

0 

AerosolModel 

0 

ModelEnvironment 

6 

TransPerKM 

0.9 

Turbulenceindex 

lE-14 

AtmSmoke 

1 

AtmSensorAltitudeIn 

1 

AtmAlpha 

1 

AtmConcLen 

1 

NumSpectDetPtsIn 

9 

Sp ect Det Wavelengt hin 1 

3.3 

Sp ect Det Wavelengt hln2 

3.5 

Sp ect Det Wavelengt hln3 

3.7 

SpectDetWavelengthInd 

4.1 

SpectDetWavelengthInS 

4.2 

SpectDetWavelengthInb 

4.5 

SpectDetWavelengthIn? 

4.6 

SpectDetWavelengthInS 

4.8 

SpectDetWavelengthIn9 

5 

NormDStarInl 

0.66 

NormDStarIn2 

0.7 

NormDStarInS 

0.74 

NormDStarInd 

0.82 

NormDStarInS 

0.84 

NormDStarInb 

0.9 

NormDStarIn? 

0.92 

NormDStarInS 

0.96 

NormDStarIn9 

1 


Units 


Comments 


0, RMS 



cm 


cm 


cm 


cm 


ft-Lamberts 


ft-Lamberts 


K 


cm 


cm 



km 


m^/g 


m^/g 


0, Beer’s Law 


Not applicable, Beer’s Law 


Not applicable, Beer s Law 



not applicable, no smoke 


not applicable, no smoke 




Same 


Same 


Same 


Same 


Same 


Same 


Same 


Same 


Same 


Same 


Same 


Same 


Same 


Same 


Same 


Same 


Same 


Same 


as FLIR92 


as FLIR92 


as FLIR92 


as FLIR92 


as FLIR92 


as FLIR92 


as FLIR92 


as FLIR92 


as FLIR92 


as FLIR92 


as FLIR92 


as FLIR92 


as FLIR92 


as FLIR92 


as FLIR92 


as FLIR92 


as FLIR92 


as FLIR92 
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Table XIX. Input Parameters for FLIR92 Part I. 


Variable 

Value 

Units 

Std Dev Optics Blur H Diam 

0.01 

mrad 

Std Dev Optics Blur V Diam 

0.01 

mrad 

H Cell Dimension 

25 

fim 

V Cell Dimension 

25 

/im 

Fields Per Frame 

2 

- 

Detectors Disp H Dir 

512 

- 

Detectors Disp V Dir 

512 

- 

Sample and Hold Filter Dir 

H 

- 

H Scene Sample Phase Ang 

0 

rad 

V Scene Sample Phase Ang 

0 

rad 

H Filter Cutoff Freq 

0 

cycles/mrad 

V Filter Cutoff Freq 

0 

cycles/mrad 

CCD H Charge Transfer Eff 

1 

- 

CCD V Charge Transfer Eff 

1 

- 

Active TV Lines Across Disp Ht 

512 

- 

Display Frame Rate 

30 

Hz 

Std Dev H Rand Image Motion Amp 

0.02 

mrad 

Std Dev V Rand Image Motion Amp 

0.02 

mrad 

Std Dev H Sinsdl Image Motion Amp 

0 

mrad 

Std Dev V Sinsdl Image Motion Amp 

0 

mrad 

Eye Limiting MTF 

1 

- 

Thrshld SNR to Resolve 4-bar Tgt 

2.5 

- 

Enviro Background Temp 

290 

K 

Det 1/f knee freq 

3 

f 

Band Avg Optics Transmission 

1 

- 

Detect H Active Dim 

23 

fim 

Detect H Active Dim 

23 

/im 

Detector Integration Time 

0.007 

s 

Tot Detectors H Direction 

512 

- 

Tot Detectors V Direction 

512 

~ 

H FIR Digital Filter Type 

none 

- 

V FIR Digital Filter Type 

none 

- 

Display Height 

0.1524 

m 

Display Luminance 

10 

ft-lamberts 

Display Viewing Distance 

0.7 

m 

Peak Det Detectivity 300K 

1 X 10^2 

cmHz2 / W 

Spectral Band Cuton A 

3.3 

/im 

Spectral Band Cutoff A 

5 

/im 













































































































Table XX. Input Parameters for FLIR92 Part II. 


Variable 

Value 

Units 

Diffraction A 

4.15 

yum 

H Focal Length 

457 

mm 

V Focal Length 

457 

mm 

Optical Aperture H Diam 

152 

mm 

Optical Aperture V Diam 

152 

mm 

FPA Readout Frame Rate 

30 

Hz 

Pathlength (sensor alt) 

1000 

m 
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APPENDIX E. ACRONYMS USED 


Table XXL Acronyms Used in this Study 


Acronym 

Name 

CCD 

Charge Coupled Device 

CTF 

Contrast Threshold Function 

FLIR92 

Forward Looking Infrared 92 

FOV 

Field of View 

FPA 

Focal Plane Array 

IFOV 

Instantaneous Field of View 

IR 

Infrared 

MDTD 

Minimum Detectable Temperature Difference 

MRTD 

Minimum Resolveable Temperature Difference 

MTF 

Modulation Transfer Function 

NETD 

Noise Equivalent Temperature Difference 

NVThermlP 

Night Vision Thermal and Image Processing 

OTF 

Optical Transfer Function 

PSF 

Point Spread Function 

SPRITE 

Signal Processing in the Element 

TAWS 

Target Acquisition Weapons Software 

TTP 

Targeting Task Performance 


101 










































THIS PAGE INTENTIONALLY LEFT BLANK 


102 



LIST OF REFERENCES 


[1] J. S. Accetta and D. L. Shumaker: editors. Infrared and Electro-Optical Sys¬ 
tems Handbook. ERIM/SPIE Optical Engineering Press, Ann Arbor/Bellingham, 
1993. 

[2] L. M. Biberman. Perception of Displayed Information. Plenum Press, New York, 
1973. 

[3] G. D. Boreman and A. E. Plogstedt. Spatial Frequency Performance of SPRITE 
Detectors. Proc. SPIE, 1309:158-166, 1990. 

[4] K. L. Davidson, T. M. Houlihan, C. W. Fairall, and G. E. Schacher. Obser¬ 
vation of the Temperature Structure Function Parameter, Cf, over the Ocean. 
Boundary-Layer MeteoroL, 15:507-523, 1978. 

[5] K. L. Davidson, G. E. Schacher, G. W. Fairall, and A. K. Goroch. Verihcation of 
the Bulk Method for Galculating Overwater Optical Turbulence. Appl. Optics, 
20:2919-2924, 1981. 

[6] R. G. Driggers, P. Gox, and T. Edwards. Introduction to Infrared and Electro- 
Optical Systems. Artech House, Boston, 1999. 

[7] G. W. Fairall, K. L. Davidson, and G. E. Schacher. Meteorological Models 
for Optical Properties in the Marine Atmospheric Boundary Layer. Opt. Eng., 
21:847-857, 1982. 

[8] R. L. Fante. Electromagnetic Beam Propagation in Turbulent Media: And Up¬ 
date. Proc. lEE, 65:1424-1443, 1980. 

[9] D. L. Fried. Optical Resolution Through a Randomly Inhomogeneous Medium 
for Very Long and Very Short Exposures. J. Opt. Soc. Am., 56:1372-1379, 1966. 

[10] J. F. Goodman. Introduction to Pourier Optics. McGraw-Hill, Inc., San Fran¬ 
cisco, 1968. 

[11] J. F. Goodman. Statistical Optics. John Wiley and Sons, Inc., New York, 1985. 

[12] A. K. Goroch. Erequency of Occurrence of Atmospheric Turbulence Conditions 
Affecting USMC Modular Universal Laser Equipment. Naval Environmental Pre¬ 
diction Research Facility, Monterey, 1980. 

[13] F. Hall. Index of Refraction Structure Parameter in the Real Atmosphere. Digest 
Printed by Optical Society of America, Boulder, 1977. 


103 



[14] U.S. Army NVESD. FLIR92 Thermal Imaging Systems Performanee Model: 
User’s Guide. Visionics and Image Signal Processing Division, Fort Belvoir, VA, 
1993. 

[15] U.S. Army NVESD. Night Vision Thermal and Image Processing Performance 
Model: User’s Manual. Modeling and Simulation Division, Fort Belvoir, VA, 
2005. 

[16] I. Overington. Vision and Acguisition. Crane, Russak and Co., Inc., New York, 
1976. 

[17] L. Scott, J. A. D’Agostino, and C. Webb. Application of 3-D Noise to MRTD 
Prediction. Proc. SPIE, 2075, 1993. 

[18] V. I. Tatarski. Wave Propagation in a Turbulent Medium. McGraw-Hill Book 
Co, New York, 1961. 

[19] R. H. Vollmerhausen and E. Jacobs. The Targeting Task Performance Metric. 
NVESD Modeling and Simulation Division, Fort Belvoir, VA, 2004. 

[20] J. C. Wyngaard. On Surface-Layer Turbulence. Workshop on Micrometeorology, 
1:101-149, 1973. 

[21] J. C. Wyngaard, Y. Izumi, and S. A. Collins. Behavior of the Refractive-Index- 
Structure Parameter near the Ground. J. Opt. Soc. Am., 61:1646-1650, 1971. 


104 



INITIAL DISTRIBUTION LIST 


1. Defense Technical Information Center 
Ft. Belvoir, Virginia 

2. Dudley Knox Library 
Naval Postgraduate School 
Monterey, California 


105 



